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This  report  is  the  result  of  a study  to  determine 
the  applicability  of  the  product  form  solution  of  Closed 
Queueing  Network  Models  for  The  DCCsystem- 1 0 . I believe  that 
this  work  will  enable  personnel  of  the  Air  Force  Avionics 
Laboratory  to  more  effectively  utilize  this  class  of  model 
for  performance  evaluation.  I hope  that  this  study  will  lay 
a foundation  for  continued  improvement  of  the  DECsystem-10 
Closed  Queueing  Network  Model. 
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McCool  and  HLchacl  E.  Price,  software  specialists  at  the 
Avioriics  Laboratory,  for  their  dedicated  support  and 
assistance.  Finally  I would  like  to  thank  Dr.  Gary  B.  Lament 
for  his  continued  support  and  advice  throughout  this  study. 

• iiarold  E.  Saxton 
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A vjlivijtion  stiKly  was  pot' t'ofmo.1  on  Closod  Quouoitii; 
Network  Models  for  the  Df!Csyst  ein- 1 0 loo.ite.l  at  the  Air  Koroe 
Aviotiios  Laboratory  (AFAi.).  Th.e  pitrpose  of  this  stady  was 
to  provide  the  personnel  of  Ai'Al.  with  an  nnderst  and  i ttjt  of 
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MODELS  FOR  THE  DECSY3TEM-10 


Chaptor  I 
Introduction 


Computer  systems  cotitinuo  to  exert  an  ever 
increasing;  influence  over  all  aspects  of  management  in  the 
Air  Force.  Over  the  last  two  deciles,  computer  structures 
and  operation  have  changed  dj-amatical  1 y . Tlie  increased 
capability  of  modern  computer  systems  has  often  been 
obtained  by  tiie  introduction  of  more  complex  systems.  Tne 
increased  computer  system  complexity  makes  it  difficult  co 
underst.ind  which  factors  most  directly  affect  th.e  operation 
of  the  system.  Tne  process  of  computer  performance 
evaluation  seeks  to  identify  these  factors. 

The  first  step  of  a computei'  per formance  evaluation 
(CFE)  is  to  identify  t!ie  measures  of  performance  wliich  are 
most  important.  These  measures  may  include  response  time, 
throughput,  resource  utilisation,  cost.  and  many  otliers 
which  are  general  in  nature  or  meaningful  only  for  specific 
modes  of  computer  operation.  on  the  desired  performatice 
measures  have  been  identified,  the  next  step  is  to  select  a 
technique  which  can  be  applied  to  the  desired  system. 

One  approach  to  computer  performance  evaluation  is 
the  use  of  modeling. 
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Ill  order  ti)  better  understand  the  v .tl  id  .it  ion  effort, 
.'1  elosed  form  .solution  of  Clo.sed  Quoueiii!^  Network  Moiel.s  l.s 
presented  with  emphasl.s  pl.ioed  on  the  physleal 
i n ter  (vot  .1 1 ion  of  m.iny  of  ttie  a.s.iumpt  ion.s  of  the  solution. 
Tlie  validation  effort.  it.selt'  eovers  multiple  worklo.ids  for 
e.ioh  eon  f ip, ur.it  ion  and  .ilso  oonslders  ,i  1 justinonts  for  the 
solution  to  compen.sate  for  as.sunpt  i on.s  which  are  not  totally 
s.itisfied  in  the  re.il  world  .system. 
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In  order  to  v.ilid.it.e  closed  qneueinp  network  models, 
.1  ihoroui'.h  knowl ..'dt',e  .)f  l.he  comput.  r systoai  .iiui  assumt'tions 
of  the  model  is  dcsir.ible.  llie  t'ollowini’,  proce. lures  were 
performed  .iuriiit;  this  stuiiy: 

I.  Tlu'  Pi'es  ys  teia- ! 0 w.i;;  in.ily.'cd  to  determine  the 

m.ijv)!'  compon.'iits  .ind  to  rel.ite  the  operat.  ion  of  l iic'ie 
c.aaponents  to  the  el.iss  of  model  t'einp,  i n v es  t.  i ,',,i  tivi  . In 
.uldition,  the  .iv.iil.ible  perform. nice  measurinp,  tools  were 
i n V e.it  i {^.it  ed  to  determine  tiie  n.itiire  of  the  d.it.i  wliicli  could 
be  use.!  t.o  v.ilid.ite  ,i  performance  iiuid.'!, 

. Tile  development  of  the  solution  for  closed 
.liieueinc,  lU'twork  models  w.i.'.  analy.’c.i  to  det.crmine  the 
r t'str  ie  t ion  .s  whioli  would  be  re.purevt  for  direi't  .ippl  i o ,it  ion 
of  this  cl.'iss  of  modi'l  t.)  t lu'  liKCsysti'm- 1 0 . 


3.  Modifications  were  made  to  the  data  structure  of 
the  Fortran  program  which  implemented  the  model.  This 
permitted  the  program  to  be  run  from  the  time-sharing 
terminals  at  the  Air  Force  Institute  of  Technology. 

M.  Several  forms  of  Closed  Queueing  Network  Models 
were  investigated  for  application  to  the  DECsystem- 1 0 . The 
results  were  summarized  in  tabular  form  and  an  analysis 
performed  on  the  observed  errors  of  each  model. 

1.4  Order  of  Pr^esentation 

Chapter  II  of  this  report  describes  the  hardware  and 
software  characteristics  of  the  DECsystem-10  at  the  Avionics 
Laboratory.  The  hardware  description  includes  the  physical 
configuration  and  the  parameters  of  tlie  individual 
resources.  Software  characteristics  include  a description  of 
the  Operating  System  resource  allocation  procedures  and  a 
review  of  the  performance  evaluation  tools  which  are 
available  on  the  system. 

Chapter  III  presents  the  development  of  a closed 
form  solution  for  closed  queueing  network  models.  The  class 
of  models  to  be  used  ia  this  report  allowr.  different  classes 
of  customers  and  four  types  of  service  center  to  represent 
the  computer  resources.  The  implications  of  assumptions 
required  for  solution  are  analyzed  as  they  relate  to  the 
DECsystem- 1 0 . 

Chapter  IV  presents  the  computationally  efficient 


5 


algorithms  for  porformatice  measures.  The  limitations  on  the 
application  of  these  algorithms  is  also  presented  in  the 
chapter . 

Chapter  V presents  the  results  of  application  of 
various  structures  of  closed  queueing  network  models  to  the 
DECsystem-1 0.  Each  model  corresponds  to  a particular 
representation  of  the  system.  The  results  of  each 
configuration  of  the  model  are  discussed  as  well  as 
comparisons  between  representations. 

Chapter  VI  presents  the  conclusions  and 
recommendations  for  future  study. 
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Chapter  H 

DECsys tem- 1 0 Computer  System 


Quantitative  measures  of  computer  system  performance 
are  only  meaningful  when  system  configuration  is  specified. 
The  system  configuration  includes  all  hardware  and  software 
characteristics.  In  general,  the  number  of  possible 
configurations  is  very  large,  even  though  the  performance 
effects  are  usually  due  to  a small  fraction  of  the  total 
number  of  parameters  (Ref  5;'-!).  The  computer  performance 
evaluation  work  by  McKenzie  (Ref  1)  includes  a 

description  of  hardware  and  software  character istic.s  of  the 
DECsystem-10  located  at  the  Air  Force  Avionics  Laboratory. 
The  addition  of  512K  words  of  core  memory  is  the  only 
significant  hardware  change  since  McKenzie's  work.  In  order 
to  avoid  excessive  references  to  earlier  works,  a brief 
description  of  the  hardware  configuration  will  be  presented 
with  emphasis  placed  on  those  areas  relevant  to  this  study. 
A more  detailed  presentation  of  selected  portions  of  the 
Operating  System  and  performance  tools  will  also  be 
presented  to  establish  a more  complete  background  for 
DECsystem-10  snodel  validation  discussed  in  Chapter  V. 

1 1 . 1 D FX_s  y s t em - 1 0 Ha r d ware 

The  configuration  of  the  AFAI.  DECsystem-10  is  shown 
in  Figure  1.  The  dual  processor  system  utilizes  two  KI10 
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central  processor  units,  which  provide  processing  service  to 
timesharing,  batch,  and  real-time  users.  Components  of 
special  interest  for  this  report  are:  the  central 
processors,  the  core  memory  and  mass  storage  devices. 

1 1 . 1 . 1 Processors 

The  KI10  processors  are  joined  in  a 
pr imary/secondary  configuration.  Both  processors  have 
access  to  all  core  memory.  However,  only  the  primary 
processor  may  perform  the  norinal  I/O  functions.  The 
secondary  processor  may  only  perform  I/O  associated  with 
real-time  processing.  The  secondary  processor  does  not 
perform  any  resource  allocation.  Therefore  the  majority  of 
the  Operating  System  overhead  will  be  associated  with  the 
primary  processor.  When  a job  running  on  the  secondary 
processor  requires  I/O,  the  job  is  stopped,  returned  to  the 
queue,  and  idontitied  for  execution  by  only  the  primary 
processor  . 

IT.  1.2 

Since  the  study  by  McKenzie  (Ref  1),  t!ie  core 
storage  has  been  expanded  with  the  addition  of  512K  words. 
The  main  memory  of  tlie  AFAI.  DECsystem-10  now  consists  of  two 
types  of  memory  units;  four  6'1K  units  of  digital  Equipment 
Corporation  (PEC)  mei.’.ory  (Ref  2 ; 70C-334-0  1 ) and  one  512K 
unit  of  Ampex  Corporation  memory  (Ref  3; 3)*  The  units 
have  the  following  parameters: 
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PM 

•Vord  Length  36  bits 

Memory  Size  5 1 2K  words 

Cycle  Time  1000  nsec 

The  addition  of  faster  memory  causes  the  required  processing 

time  for  a job  to  decrease.  The  amount  of  decrease  is 
dependent  upon  the  structure  of  the  job  and  the  number  of 
memory  references. 

II . 1 . 3 2^931^8? 

Two  disk  systems,  RP03  and  RP04,  are  used  by  the 
DLCsystem-IO.  Four  RP03  disk  units  are  interfaced  through 
the  MX10  memory  multiplexor  channel,  while  four  high  speed 
RP04  disk  units  are  connected  thr’ough  an  RH-10  control  unit. 
These  disk  units  have  the  following  characteristics  (Ref 
2:70C-384-01g) : 


_RP03 

RP0_4 

number  of 

Units 

4 

4 

Capacity 

(Million  Words) 

10.24 

20.00 

Transfer 

Rate  (words/ sec) 

66 , doY 

178,971 

Average  Access  Time  (sec) 

4 7 . 5 

36.8 

Of  special  concern  to  this  study  is  thi?  system  use  of  each 
disk  type.  One  RP03  unit  is  used  for  system  files,  while  one 
RP04  unit  is  used  for  swapping.  Under  heavy  load  conditions, 
interference  may  result  for  either  disk  type. 


AMPEX 

36  bits 
512K  words 
920  nsec 
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II.1.U  Real-Time  Pfooer.r. inc  Krf<^ols  oti  ijaixlvvare 
Conf  if'urat  ij3ri 

McKenzie's  work  provides  a compreliensi ve  discussion 
of  pr'oblems  arising;  from  real-time  programming  applications. 
The  excessive  swapping  overhead  induced  by  real-time 
program's  use  of  memory  has  been  nearly  eliminated  by  the 
additional  core  storage.  ilowever,  future  applications  may 
rejuire  increasingly  large  amounts  of  core  and  again 
restrict  time-sharing  user's  access  to  memory.  Another 
consequence  of  real-time  processing  is  the  allocation  of  a 
single  processor  to  time-sharing  users.  The  dedicated  use  of 
the  secondary  processor  is  investigated  further  in  this 
study . 

1 1. 2 DKCsyste:n-  1 0 Operatuig  System 

The  TOPS-10  operating  system,  or  tlOtJlTOK,  performs 
the  accounting,  scheduling,  resource  allocation,  and  service 
routines  necessary  to  operate  in  a inu  1 1 i pr ograiiiming  , time- 
sharing environment  (Ref  't ) . 'Uimerous  tables  and  queues 
are  maintained  by  the  system  in  order  to  perform  the 
required  functions.  The  MONITOR  maintains  control  of  all 
processing  by  responding  to  a clock  Interrupt  generated  for 
e.ich  processor.  The  clocks  are  run  at  line  frequency.  The 
interrupts  therefore  bracket  a time  slice  1/bO  second  in 
length.  This  measure  of  time,  called  a jiffy,  will  pi'ove 
useful  when  discussing  system  operation. 


1 1 


There  are  a variety  of  module.s  and  programs  within 
the  Operating  System  which  permit  it  to  perform  its  required 
tasks.  Among  these  are  the  Conti'ol  Cycle  modules,  including 
resource  allocation  modules,  and  various  evaluation  tools 
which  enable  system  performance  measurements  to  be  obtained. 


1 1. 2.1 

Beginning  with  the  clock  tick  interrupt,  every 
sixtieth  of  a second,  the  cyclic  routines  of  the  MONITOR  are 
executed  and  then  control  is  returned  to  the  user  program. 
Ttie  activity  of  the  control  cycle  takes  on  added  importance 
for  this  study.  One  of  the  primary  data  gathering  programs, 
METER,  utilizes  code  located  in  strategic  tnodules  within  the 
MONITOR.  The  major  MONITOR  modules  activated  during  each 
clock  tick  are;  CL0JK1,  COMCON  and  SCHEBI.  The  sequence  of 
these  modules  is  shown  in  Figure  2. 


Figure  2.  Control  Cycle  Sequence 
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II.2.  1 . 1 


Clocki 


r 

I 

i 

I 

I 

I 

I In  CLOCKI,  the  basic  accounting  functions  of  the 

* MONITOR  are  performed.  If  a user  program  was  running,  the 

j accumulated  time  since  the  last  accounting  update  is  added 

: to  the  job's  run  total  in  the  MONITOR'S  Process  Data  Block. 

' When  there  are  no  Jobs  ready  to  run,  a Monitor  routine 

I 

: called  the  Null  program  is  given  control  of  the  processor. 

If  the  Null  program  was  running,  the  time  is  added  to  the 
total  Null  job  time  in  the  CPU  Data  Block.  If  the  Null 

B program  was  running  and  there  was  at  least  one  job  in  any  of 

3 

1/  the  run  queues  which  could  not  bo  run  for  some  reason,  then 

''  the  time  is  considered  "Lost"  CPU  time.  One  situation  where 

Lost  time  may  occur  is  when  a job  has  been  swapped  out  to 

disk  and  is  scheduled  to  run  before  it  can  be  swapped  back 

into  memory.  If  no  jobs  were  in  a run  queue  while  the  Null 

I program  was  processed,  the  time  is  considered  "Idle"  time 

I 

I (Ref  ’* ) . This  time  is  used  to  calculate  the  CPU 

utilization.  In  addition  to  the  accounting  function,  CLOCKI 

I 

I is  entered  at  the  end  of  the  control  cycle  to  perform 

J context  switching  if  needed.  Context  switching  is  the 

process  of  saving  user  register  coi.t^ents  and  pointers.  This 
is  required  only  if  a new  job  has  been  selected  by  the 
scheduler.  When  the  accounting  functions  are  complete, 
CLOCKI  calls  a routine  which  determines  if  there  are  any 
commands  waiting  to  be  processed.  If  there  i.re  commands 
waiting,  control  is  passed  to  tlie  command  processoi'  rcuti..  . 
COMCON. 


13 


11.2.1.2 


Com:natKi  rfoot'S-sor 


The  Command  Processor  is  the  MONITOK  module 
responsible  for  interpret  int;  the  user's  retjuest  from  the 
terminal  and  then  passing;  control  to  the  proper  routine  to 
satisfy  the  command.  After  the  user  lias  entered  a command 
from  the  terminal,  the  command  is  stored  in  a buffer 
internal  to  the  MONITOR.  If  there  are  commands  in  the 
buffer,  then  COMCON  will  activate  routines  which  will  gather 
the  data  necessary  for  later  reference  to  the  appropriate 
user's  terminal.  The  command  processor  will  not  process  more 
than  one  comm.ind  during  each  control  cycle.  If  more  than  one 
command  is  waiting,  the  command  processor  will  be  called 
daring  successive  cycles  to  process  the  remaining  commands 
(Kef  ) . This  is  consistant  with  the  basic  MONITOR 
philosophy  which  is  to  run  the  control  cycle  quickly  to 
cor.iplot  iv)n  and  allow  the  user  programs  tlie  majority  of  the 
pi'ocessing  timi?.  gome  commands  may  require  more  time  than  is 
desirable  to  be  taken  during  the  control  cycle.  These 
commands  are  set  up  as  special  MONITOR  routines,  and  are  run 
during  the  user's  time  slice  as  one  oi'  the  user's  Jobs. 

II..:.  1.3  SCIU'DI 

The  scheduler,  SCHKl'l,  is  the  MONITOR  module 
responsible  foi*  controlling  the  .vl'  'cation  of  system 
resources  among  user  jobs.  The  Ollvisy.s  i 10  is  a 
multiprogramming  system  which  permits  numerous  vis.  • ' to 
reside  in  core  with  concurrent  operation.  It  is  the  functio 


of  the  scheduler  to  determine  which  job  will  run  during  each 
time  slice.  When  swapping  occurs,  as  is  permitted  at  AFAL, 
then  jobs  may  reside  in  core  or  on  secondary  storage  devices 
such  as  the  RP04  disks.  Because  of  this  situation,  the 
scheduler  must  not  only  decide  which  job  will  run,  but  also 
call  a Swapper  routine  to  determine  which  jobs  will  be 
swapped  to  or  from  secondary  storage  and  when  this  swapping 
will  occur  . 

The  overall  objectives  for  the  scheduler  are: 

1)  Provide  for  sharing  the  machine's 
resources  among  the  user  jobs  so  that  each 
receives  his  share  on  a timely  basis. 

2)  Provide  fast  response  capability  for 
jobs  that  require  conversational 
interaction. 

3)  Provide  a very  fast  response 
capability  for  jobs  performing  tasks  that 
require  "real  time"  interaction. 

4)  Make  efficient  use  of  all  system 
resources  (Ref  7 ; 1 ) . 

To  meet  these  objectives,  the  scheduler  uses  a 
modification  of  the  Round  Robin  (RR)  algorithm,  one  of  the 
most  widely  used  scheduling  algorithms  for  time-shared 
computer  systems.  The  structure  for  tl'.c  Round  Robin  system 
is  shown  in  Figure  3-  Newly  arriving  customers  take  their 
place  at  the  back  of  the  queue  and  progress  to  the  head  of 
the  queue  in  f irst-come-first-serve  fashion  (Ref  8; 166). 
When  a job  is  ready  to  receive  service,  it  is  allocated  a 
service  limit  called  a quantum.  If  the  job  still  requi  ' s 
service  when  the  quantum  has  expired,  it  is  returned  to  the 
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opportunity  to  use  the  rosouroe  ansociotod  with  the  queue. 

Kaeh  of  the  jobs  receive  a "fair  share"  of  the  resource 
capacity.  In  this  instance,  "fair  share"  means  that  no  job, 
rct^ardless  of  its  service  requirements,  may  take  over  the 
resource.  Diffei'ent  users  may  have  ilifforinj;  opinions  on  j 

what  is  a "fair  share"  for  their  particular  jobs.  Some  jobs 
may  require  mucli  less  service  than  would  be  allocated  to 
them  in  the  Hound-I\obin  alpjOrithm.  for  these  jobs  it  might 
be  considered  unfair  to  require  them  to  wait  in  the  queue 
for  large  time  quantum  which  they  will  not  fully  utilize. 

Interactive  jobs  are  one  class  of  job  which 
generally  have  small  processing  retiui  rements . Tlie  Text 
Kditor  and  Corrector  (I’KCO)  is  an  example  of  an  interactive 
job.  For  this  type  job,  a fast  response  time  capability  i-s 
desircvi.  when  the  user  enters  a TKCO  commanvi  , he  expects  a 
quick  response,  ([)roforably  less  than  5 seconds).  The  time 
required  to  cycle  through  the  I\).?  queue  will  be  highly 
dependent  upon  the  number  of  jobs  in  the  queue.  If  the 
system  is  heavily  loaded,  the  time  required  to  cycle  through 
P,).?  could  easily  be  several  times  the  desii'od  I'csponso. 

To  meet  this  requirement  for  faster  response  time 
for  small  jobs,  another  queue,  PQ1,  is  UM.d.  T!iis  iiuoue  will 
have  a higlier  priority  tlian  PQ2.  A time  iiuanium  is  also 

i 

established  for  jobs  in  Pi.)  I . However,  tae  FO I quantum  will 
t'enerally  be  much  smaller  than  tlKit  received  in  PJ? 

(Kef  7;2).  If  a job  requires  mare  service  aftc'r  its  PQl 
quantum  has  expired,  then  it  will  b('  requeued  io  the  back  of 
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the  PQ2  queue.  In  this  manner,  all  jobs  receive  a small 
burst  of  service  while  in  the  PQ1  queue.  Jobs  which  still 
require  service  enter  the  Round  Robin  cycle  for  PQ2,  where 
they  remain  until  completed. 

Certain  programs  may  require  even  faster  response 
and  better  performance  than  is  obtained  by  using  PQ1.  To 
satisfy  this  class  of  programs,  additional  high  priority 
queues  are  established.  They  are  labeled  liPQI  through  liPyi5 
(Ref  7 ; 6 ) . The  Line  Printer  Spooler  and  certain  computer 
performance  evaluation  programs  are  examples  of  the  type  of 
job  which  require  these  queues.  If  not  in  use,  the  Line 
Printer  Spooler  is  swapped  to  secondary  storage.  When 
pr’inting  is  required,  the  spooler  must  be  swapped  into  core 
ahead  of  other  user  jobs  and  must  have  access  to  the  CPU 
anead  of  user  jobs  so  that  the  output  buffer  may  be  filled. 
Performance  measuring  programs  often  require  rapid  access  to 
MONITOR  data  tables.  If  liigh  priority  queues  were  not 
available,  the  CPL  programs  would  have  to  wait  in  the  PQ1 
queue  and  possibly  lose  data.  Jobs  in  the  High  Priority 
queues  may  acquire  all  system  resources  ahead  of  user 
programs.  However,  tliose  jobs  which  are  permitted  into  HPQ 
are  generally  1/0  bound  jobs  so  that  most  of  the  CPU 
capacity  is  available  to  the  user  jobs. 

The  final  objective  of  the  scheduler  is  to  operate 
the  system  as  efficiently  as  possible  within  the  previously 
itientioned  performance  constraints.  One  approach  is  to 
balance  the  CPU  and  I/O  bound  jobs  in  core  so  that 
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multiprogramming  is  made  most  efficient  by  ■ overlapping  I/O 
with  processing.  Wlien  the  primary/secondary,  dual  processor 
configuration  is  used,  additional  considerations  must  be 
included  in  the  scheduling  philosophy.  Since  the  secondary 
processor,  CPU1,  has  no  I/O  capability  for  time-sharing 
users,  it  would  be  desirable  to  process  CPU  b.ound  jobs  on 
that  processor.  That  approach  would  permit  I/O  bound  jobs  or 
jobs  with  small  processing  requirement.s  to  be  run  on  the 
primary  processor.  This  approacli  is  incorporated  into  the 
scheduler  so  that  consideration  is  given  the  type  of  job 
assigtied  to  each  processor.  Jobs  to  be  run  on  CPUO  are 
selected  in  order  from  the  front  of  the  higest  priority  run 
queues  as  follows:  1IPQ1  thi'ough  UPQ15,  PQ1,  and  finally  PQ2. 
Using  the  assumption  that  jobs  which  have  already  run  for  a 
relatively  long  period  of  time  will  be  CPU  bound,  the  jobs 
to  be  run  on  CPU1  are  selected  from  tlie  front  of  the  lower 
priority  run  queues  in  the  following  order:  PQ2,  PQ1,  and 
finally  the  HPQ's.  Another  factor  is  the  realization  that  a 
compromise  must  be  achieved  between  throughput  and  short 
term  response.  One  of  the  controlling  considerations  for 
this  compromise  is  the  number  of  jobs  in  the  short  time  run 
queue,  PQ1,  and  the  long  run  queue,  P02.  The  number  of  jobs 
in  each  queue  is  directly  dependent  upon  the  size  of  the 
quantum  allocation  in  each  queue. 

11.2.2  Quantum  Run  Time 

The  flow  of  jobs  through  the  processor  run  (pieues  is 
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greatly  affected  by  the  amount  of  service  they  receive 

before  being  requeued  and  also  by  the  length  of  time  they 
are  permittv  to  remain  in  main  memory.  The  job's  quantum 
run  time  is  decremented  while  the  job  is  being  processed. 
The  in-core  protect  time  is  decremented  if  the  job  is 
scanned  by  the  scheduler,  whether  or  not  it  is  being 

serviced.  The  in-core  protect  time  affects  processing  of 
jobs  in  another  manner  as  well.  When  the  in-core  protect 
time  expires,  the  job  is  requeued  to  the  end  of  PQ2.  This 
may  occur  well  before  its  alloted  time  quantum  has  expired, 

and  even  if  the  job  is  not  actually  swapped  to  secondary 

storage . 

Quantum  run  time  is  computed  by  the  following 

formula: 


QUANTUM  HUN  = M IN  [ QMX  , QAD+ ( S IZ  t;»QML) /QR  ANG  ] 

Size  is  the  job's  size  in  K (1024  words).  QMX,  QAD,  QML  and 
QRANG  are  Operating  System  parameters  which  are  stored  in 
MONITOR  tables  with  values  specified  for  each  queue  and  with 
the  option  to  change  the  values,  depending  upon  which 
swapping  device  is  used.  For  the  synthetic  jobs  used  in  this 
study,  the  quantum  run  allocations  are:  3 jiffies  for  PQ1, 
and  30  jiffies  for  PQ2. 

A similar  formula  is  used  to  determine  in-core 
protect  time.  In  this  case  the  Minimum  Core  Usage  Function, 
MCUF,  is  defined  as  follows: 

MCUF  = (PROT  » SIZF)  + PROTO 
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PROT  and  PROTO  are  intornal  Monitor  parameters.  The  current 
parameters  at  AP'AI.  produce  an  in-core  protect  time  value  of 
2 seconds,  120  Jiffies,  for  the  synthetic  jobs  of  this 
study . 

The  determination  of  quatitum  run  time  for  each  queue 
and  in-core  protect  time  can  have  a profound  effect  on 
system  performance.  for  P01 , the  quantum  run  time  is  a 
measure  of  the  time  the  Job  receives  rapid  service  after 
entering  tiie  system.  Wiien  this  quantum  expires,  the  Job  is 
requeued  to  the  end  of  PQ2  where  it  receives  a new  run 
quantum.  However,  the  Job's  in-core  protect  has  been 
decreased.  This  procedure  permits  ;:Mall  fast  Jobs  to  receive 
exceptional  service  scheduling  response  and  also  reduces 
swapping  since  the  longer  Jobs  are  still  permitted  to 
receive  some  service  in  PiJ2  before  they  are  swapped  out 
(Ref  4) . 

For  Jobs  in  P'J2,  the  quantum  run  time  and  in-core 
protect  time  control  the  bias  of  the  scheduler  for  I/O 
versus  CPU  Jobs.  If  both  parameters  are  increased,  there 
will  be  a decrease  in  thi?  number  of  Jobs  whoso  time  slice 
has  expired.  Accordingly,  swapping  decrea.se  and 
throughput  should  improve  (because  of  reduced  overhead  in 
swapping).  Respoiise  is  degraded  in  this  instance  because 
Jobs  will  tend  to  wait  longer  to  swap  into  core.  Decr'casitig 
the  two  parameters  will  h.ave  the  opposite  effect  on 
throughput  and  response.  The  ratio  of  the  two  parameters  is 
also  important.  Loni;er  cpiantum  run  times  favor  CPU  Jobs, 
while  longer  in-core  protect  times  favor  I/O  Jobs. 
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Two  odd  itioiiol  rostr  io  t ions  afo  plaood  on  swapping; 

(Ref  4): 

1)  Real-time  jobs  may  not  be  swapped 
out  of  core. 

2)  Jobs  with  incomplete  I/O  request  may 
not  be  swapped  out  of  core. 

Tlie  scheduler  calls  up  the  swappint;  routine  which  selects 
jobs  to  be  swapped  into  or  out  of  core,  within  the 
restrictions  stated  above.  Jobs  to  be  swapped  in  are 
selected  from  t!ie  front  of  the  run  queues  in  order  of 
priority:  liPQ's,  PO 1 , and  POP.  Jobs  to  be  swapped  out  are 
selected  from  the  end  of  the  run  queues  from  lowest  to 
hi^jhest  priority:  POP,  PO 1 , and  llPO's. 

1 1.2. 3 System  Tools 

Tliree  performance  measuring,  prot’rams  .ire  available 
on  the  DKCsystem- 1 0.  They  are  SYSi'AT,  TRACK,  and  MKlbiR.  In 
addition,  the  SCRIPT  proi^ram  per:iiits  the  simultaneous 
execution  of  several  copies  of  a sitp’le  proi;:'am.  execution 
of  various  combinations  of  pi'op.rams  .ilonp,  witti  one  of  the 
performance  measurini.;  prot;rams  allows  the  collection  of 
detailed  information  on  system  per  foi'm.ino  e under  controlled 
cond  i tions  . 

1 1. 2. 3.1  SCRIPT 

SCRIPT  is  a system  pro{:,r.im  which  fiermits  numerous 
copies  of  a progr.im  to  be  executed.  For  tliis  study,  two 
basic  programs  were  used  as  syntlietic  workload.  Using 
SCRIPT,  up  to  14  copies  of  each  rror.r.im  can  be  logged  ont.c 
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psoiKio  toloLypos.  Iti  >uldLtion,  one  ol'  the  pee  t'oniKHice 
ne.isui'iiig  proi^cmns  can  be  initiated  each  time  so  that  the 
resulting  performance  data  would  be  a measure  for  a known 
workload.  Simultaneous  execution  of  several  "scripts"  permit 
viny  combination  of  jobs;  li:nited  only  by  the  system's 
c apac i ty  . 

I I.  2.  1.2  Systat 

Tlie  SYST.'ir  program  provides  the  user  or  computer 
eeiter  manager  with  getu'ral  system  status  information. 
Befj;\  loj’t’ing  ctue  the  systei.i,  a user  may  use  the  SYSTAT 
command  to  di'terpiine  tlie  number  of  jobs  on  the  system  and 
t hi'  pereeiit  Null  time  (idle  time  plus  lost  time).  Limited 
I n f jfiiat  ion  a'ooui  ea.  h.  usei*  is  also  available.  The  SYSTAT 
pf'/t  :n  nrovid.es  oao’!  user's  job  number,  [irogram-projecL 
nu:;ibet  , and  I '.ji.!  ’t,\'..ee.  Also  available  is  the  name  and 

s i ae  of  tlie  pro/.ra.u  each  user  is  I'unninj'  and  the  aeeumi;  I at  ed 
!'un  rime  since  lu.'f^ing  ontv>  td'.e  s\stei'i.  Another  ineiisurc  of 
systc'i  load  whieli  is  provivled  by  SYSlAT  is  a listing  of 
active  devices  and  amount  of  fre<.'  slorage  for  eacli  disk 
structure.  Of  specr-il  interest  to  the  system  manager  is  the 
a,.tount  of  swapping  space  used,  virtual  core  used.,  swapping 
ratio,  virtual  core  saved  bv  sharing,  and  averag.e  job  .size 
(licf  10).  S'Ybi'Al'  does  provide  useful  i n foi'm.tl  i on  under 
actual  lo.id  eondili  'ns.  However,  a more  flexible  program  is 
needed  to  gather  tlte  .lat..i  nee. led  for  model  v .j1  i d at- i on  . Tlie 
desired  flexibility  is  provided  by  the 
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TilAOK  progr.im. 
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T r ao  k 
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The  TRACK  procrain  may  be  used  to  monitor  the 
progress  and  performance  of  individual  jobs  or  it  may  be 
used  to  measure  the  performance  of  the  entire  system.  For 

this  study,  the  TRACK  program  was  used  only  to  gather  data  | 

on  system  performance.  The  program  performs  a "peek"  at  j 

selected  performance  pv^rameters  at  intervals  specified  by 

ttie  user.  TRACK  is  initiated  by  entering  a command  string 

which  specifics  requirements  such  as:  typo  of  report,  use  of 

high  priority  queue,  look  in  core,  length  of  peek  interval, 

and  number  of  intervals  used  in  each  report  (Ref  9). 

The  report  generated  by  TRACK  will  include  a 

frequency  distribution  for  most  parameters  measured  as  well 
as  the  average  value  and  standard  deviation  for  each 
parameter.  Information  concerning  CPU  performance  was 
especially  important  to  this  study.  The  TRACK  information 
for  each  processor  includes:  percent  Null  time,  percent  Lost 

I time,  and  percent  time  spent  for  overhead.  Additional 

; measurements  useful  for  model  validation  are  the  number  of 

I , 

Jobs  in  the  run  queues,  and  the  number  of  jobs  in  T/0  wait 

queues  other  than  teletype.  ^ 

1 1 . <? . 3 • 't  Meter 

The  meter  facility  permits  t(ic  implementation  of 
"meter  points"  throughout  the  MONIIOR,  wherever  a value  or 
event  of  potential  interest  mav  occur.  The  Meter  facility  is 
u.sed  to  select  and  collect  performance  statistics  from 

I 

I 

.?4  . 
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MOMITOl?.  It  allows  priv  Helped  usof-:nodo  pr'oj',r:ims  to  collect 
system  data  for  performance  analysis  or  tunitiii  of  jobs  or  of 
the  system.  The  Meter  facility  is  capable  of  providintj 
simple  statistics  or  raw  data.  The  raw  data  outputs  require 
processing  and  formatting  after  the  program  has  been  run. 
however,  data  in  this  form  provides  a more  detailed  picture 
of  the  system  operation  and  is  more  desirable  for  the  early 
stages  of  Cri£  model  development. 

The  Meter  facility  is  comprised  of  three  mechanisms 
which  are  used  for  tlie  collection  and  disposition  of 
performance  data  (Ref  9 ; 1 . 1 ) : 

1)  Meter  Points  are  short  sections  of 
code  (often  only  two  instruc tions)  located 
throughout  the  MORTTOR.  This  code  test  for 
activation  of  its  particular  sec’'ion  of  code 
and  if  active,  calls  a Meter  Point  Routine. 

?.)  Meter  Point  Routnies  are  sliort,  fast 
routines  wiiich  process  the  meter  values 
before  passing  control  to  the  Meter  Channel. 

T)  Meter  Channels  are  fast  routines 
which  dispose  of  the  processed  values  into  a 
buffer  under  user  control. 

The  buffer  used  to  aoeumulate  metered  data  is  located  in  the 
metering  job's  core  image.  Tlie  meter  i nr,  job  uses  the 
hibernate  monitor  call  to  synchronize  t!ie  acouinul  ation  and 
processing  of  data.  The  buffer  length  is  set  to  an  integral 

J 

power  of  two  so  that  indexing  may  be  performed  by  a simple 
modulo  operation.  The  MOIJITOR  does  not  guard  against  buffer 

overflow,  but  does  provide  the  hibernate  function,  kdien  ^ 

large  amounts  of  data  are  required,  the  motoring  job  may 
have  to  be  placed  in  a lil;h  priority  queue  in  order  to 

prevent  buffer  writeover.  i 

L.  ...  i 


The  data  gatherini’  routines  used  in  tliis  validation 


study  included  three  meter  points: 

1)  Meter  point  1,  Job  State  Queueing. 
This  point  is  located  in  SCHKU1  at  the  point 
where  destination  queue  has  bee;;  determined. 

2)  Meter  point  2,  Context  Switching. 
This  point  is  located  in  CL0CK1  where  it  has 
determined  that  a different  job  is  to  be 
run  . 


3)  Meter  point  500,  Wait  State  Code. 

This  point  is  located  in  SCI!K01  where  a 
change  in  state  code  has  been  required. 

Wlien  raw  data  is  collected,  it  must  be  received  and  buffered 

as  it  is  generated.  In  order  to  avoid  losing  data,  a routine 

which  ' locked  in  memory  must  be  provided  so  that  tlie  data 

may  be  dumped  onto  secondary  storage.  If  the  routine  is  not 

locked  in  memory,  it  may  not  be  swapped  back  into  core  in 

time  to  make  the  transfer  when  the  system  is  experiencing  a 

high  swapping  rate. 


1 1 . 3 Sunim rir  y 

The  DECsystem-IO  is  a dual  processor  computer  system 
with  a complex  Operating  System  which  permits  time-sharing, 
batch  and  real-time  users  to  share  the  resources  of  the 
system.  The  complexity  of  the  system  ari.ses  not  so  much  from 
the  hardware  configuration  as  from  the  Operating  System 
which  attempts  to  both  optimize  resource  performance  and 
provide  so::ie  degree  of  equality  to  the  users. 

The  system  measuring  tools  provide  various  levels  of 
information  which  may  be  useful  for  performance  evaluation. 


i 
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A desirable  goal  of  a performance  model  is  that  it  require 
only  readily  available  system  information  in  order  to  make 
acceptable  predictions  of  the  system  performance.  In  order 
to  fully  develop  the  model,  it  is  often  necessary  to  obtain 
detailed  information  of  the  system.  The  DECsystem-10 
provides  evaluation  tools  whicli  present  both  levels  of 
information  which  greatly  aided  this  validation  effort. 


Chap'-^  I 


Closed  Queue  ing  Network  M_od^^l 

The  paper  by  McKenzie  (Ref  1)  used  several  Closed 
Queueing  Network  Models  to  evaluate  the  performance  of  the 
DECsystem-10  at  the  Air  Force  Avionics  Laboratory  (AFAI.).  In 
order  to  validate  models  of  this  type  under  various  workload 
conditions,  a thorough  understanding  of  the  assumptions  and 
solution  techniques  of  the  model  is  desirable. 

The  solution  of  the  Closed  Queueing  Network  Models 
used  in  this  report  is  based  on  the  work  of  Baskett,  Chandy, 
Muntz,  and  Palacios  (Ref  11).  Tliis  class  of  model  is  an 
extension  of  earlier  works  by  Jackson  (Ref  12)  and  Gordon 
and  Newell  (Ref  13).  The  solution  technique  of  Gordon  and 
Newell  is  presented  in  Appendix  A.  In  this  chapter  the  work 
of  Baskett,  et.  al.  is  reviewed  to  provide  the  background 
for  a discussion  of  the  assumptions  which  most  directly 
affect  the  application  of  tliis  class  of  model  to  the 
DECsystem- 1 0 . 

1 1 1 . 1 Work  I^a  sk e 1 1 

Baskett,  et.  al.  have  extended  earlier  works  to 
allow  different  classes  of  customers.  A represcntatioti  of  a 
Closed  Queueing  Network  Model  is  shown  in  Figure  'I.  The 
introduction  of  multiple  customer  classes  overcomes  the 
limitation  of  previous  works  whio'.i  required  that  all 
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c u.':  Loiiit'fii  belong  to  the  n.itiio  olns.s.  Cu  sLoiium's  feceivo 

sofv  Loo  ;il  the  rosourcos  in  tlio  nolwoi'k  neooc'j  Ing  to  nn 

oxponontinl  not'vico  time  il  istr  ibut  ion  specified  for  eacli 

node  in  the  network.  Me;in  service  time  for  clnss  r 

eustoiiu-rs  ;it  node  i is  given  by  I/mj,.  Movement  of  customers 

between  resources  is  specified  by  a fixed  transition  matrix 

l'=Lp  . . ] wliere  p . is  the  probability  that  a class  r 

I r ; j s 1 r ; j ;> 

eustoiier  leaving  the  Ith  ."aM’vice  center  will  change  to  a 
class  s customer  uui  g.o  to  the  jth  service  center.  Customers 
within  a given  class,  are  described  by  the  s.ime  service  time 
distribution  and  tlie  same  transition  probabilities.  The 
emphasis  of  this  pajH'r  is  towards  validal.ion  of  the  moilel 
for  custoiiu'rs  which  do  not  ciiang.e  class.  I.atei'  eominit  at  iona  I 
coils  ider  at  ion  will  expresr.ly  prohibit  class  chang.c-'.  This 

sim[>ly  means  that,  p . . when  r/ts. 

1 r ; J 

The  assumption  of  exponential  service  time 
vi  istr  ibut  ions  is  not  valiii  for  all  lu'sources  wliii'h  are 
likely  to  appear  in  a comimtcr  system  model.  Ai'tual  system 
measurements  on  the  Micliigan  id*  rm  i n.i  1 System  indicated  that 
tlie  s.erv  ice  times  for  that  system's  iiwapping,  drum  were 
ii  y pere  X ponen  t i a I (Kef  PI;PI7). 
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which  iiiay  be  used  to  represent  the  service  time  distribution 
of  a service  center.  The  probability  that  a customer 
receiving  service  at  the  ith  center  is  in  the  1th  stage  of 
service  is  given  by  (Ref  11): 

1 

"in  = n ‘‘I'-J 

J=1 

Using  the  concept  of  classes  of  customers  and  tlie 
method  of  stages,  Paskett,  et.  al.  present  the  equilibrium 
state  solution  for  a queueing  network  with  four  allowable 
types  of  service  center.  These  four  types  of  service 
centers  are  (Ref  11): 

Type  1:  Single  server,  Fi rst-Come-Fi rst-Serve  (FCFS) 
queueing  discipline,  all  classes  must  have  the  same 
exponential  service  time  distribution. 

Type  2:  Single  server,  Processor  Sharing  (PS) 
queueing  discipline,  different  classes  of  customers  may  have 
different  service  time  distributions.  All  service  time 
d istr  ibut  ions  must  liave  rational  Laplace  Transform. 

Type  3:  Infinite  servers  (IS),  no  customer  waits  for 
service  since  th.ere  are  at  least  as  n;.iny  servers  as 
customers.  Pi  f fere’'.'i..  classes  of  custodiers  may  have  unique 
service  time  distributions  so  Umii;  as  the  distributions  liave 
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rational  Laplace  Transforms. 
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i.  Type  U:  Single  server,  Preemptive-Resume  Last-Conie- 

f 

; First-Serve  (LCFS)  queueing  discipline,  each  class  may  have 

f a distinct  service  time  distribution  with  rational  Laplace 

Tr  ans  form . 

j 

Service  centers  with  multiple  servers  and  queue  dependent 
service  rates  will  be  considered  after  tlie  development  of 
the  solution  for  the  network  of  single  servers  centers. 

1 1 1 . 2 Network  SJ:.  a t_o  o n r esent a t i_o n 

The  state  of  a model  may  be  described  at  various 
levels  of  abstraction.  For  example  the  state  may  be 
specified  by  the  total  number  of  customers  in  the  system, 
the  number  of  customers  at  each  resource,  or  a detailed 
description  of  the  ordering  atui  characteristics  of  the 
customers  at  each  resource. 

Fiaskett  et.  al.  represent  the  state  of  the  model  by 
a vector  ( ,X2,...,Xk|)  where  Xj^  represents  the  condition  of 
thi  ith  service  center.  The  method  of  describing  a center's 
condition  ilepends  on  the  type  of  center. 

If  the  service  center  is  FCFS  then 
^ ,Xj^2  t • • • >^in^  ^ 

where  n^  is  the  number  of  customei’s  at  center  i and  x^  j is 
the  class  of  customer  whicFi  is  in  the  jth  position  of  the 
FCFS  order.  The  single  server  center  is  considered  at  this 
point.  TFierefore,  only  the  first  customer  in  the  queue  is 
served  while  tFie  remainder  are  waiting. 

3^ 


If  the  service  center  is  either  PS  or  IS,  the 


center's  condition  is  described  by: 

Xi= ( Vi ^ ,Vi ^ , . . . , Vip) 


where  is  a vector  (m ^ ^ • ,m^  ^) . The  Ith  component 

of  V . , is  the  number  of  class  r customers  in  center  i and 
1 r 

in  the  Ith  stage  of  service.  Each  class  of  customer  may  have 
its  service  time  distribution  represented  by  a different 
number  of  stages.  Therefore,  u^^  is  the  number  of  stages  for 
class  r customers  at  tlie  ith  center.  All  customers  in  the  PS 
center  continually  receive  some  amount  of  service  from  the 
single  server.  The  percentage  of  full  server  capacity  which 
each  customer  receives  is  dependent  upon  the  total  number  of 
customers  at  the  center.  Customers  at  an  IS  service  center 
receive  full  server  capacity  since  there  are  at  least  as 
many  servers  as  customers. 

If  the  service  center  is  governed  by  Preemptive- 
Resume  LCFS  discipline,  then  the  center's  condition  is 
described  by: 


X t r I ,m  p , ( r^  .m^)  ,...,(  r^  _ ''''n  • ) ^ 


where  n^  is  the  number  of  customers  at  the  ith  center  and 

(r.,m.)  describes  the  jth  customer  in  LCFS  order;  r.  is  tlie 

class  of  the  customer  and  m,  is  the  stage  of  service  the 

J 

customer  occupies. 

The  state  description  for  a closed  queueing  network 
with  multiple  classes  is  much  more  complex  than  that  used 
for  the  single  class  moilel  described  in  Appendix  A.  For  tlie 


For  the 


system  with  multiple  classes,  the  state  description  must 
consider  all  partitions  of  the  customers  among  the  service 
centers  and  the  ordering  and  characteristics  of  the 
customers  at  each  center.  The  concept  of  an  equilibrium 
condition  however,  is  very  similar.  The  equilibrium  state 
probabilities  must  satisfy  the  following  (Ref  11): 

V states,  P(Sj)[rate  of  flow  from  Sj  to  Sj^]  = 

all  states  P(S^)[rate  of  flow  out  of  S^] 

3 . 
j 

Chandy  (Ref  16)  calls  these  equations  the  global  balance 
equations.  He  also  defines  a set  of  local  balance  equations 
which  require  that  each  term  on  the  right-hand  side  be  equal 
to  a particular  subset  of  terms  on  the  left-hand  side.  Thus, 
"a  local  balance  equation  equates  the  rate  of  flow  into  a 
state  by  a customer  entering  a stage  of  service  to  the  flow 
out  of  that  stage  due  to  a customer  leaving  that  stage  of 
service"  (Ref  11;252).  The  examples  presented  by  Baskett, 
et.  al.  specify  a local  balance  equation  for  each  class  of 
customer  and  each  service  center  where  each  center  has  only 
one  stage  of  service. 

The  solution  for  the  multiple  class  model  yields  a 
set  of  simultaneous  equations  similar  to  those  found  for  the 
single  class  model.  These  are  presented  by  Wong  as 
(Ref  17): 
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j = 1 ii=1 


e 

ir 


e . P . 
js  js;ir 
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r = 1 , 2 . . . R 


(2) 
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v.'here  is  oquiviilent  to  in  the  notation  used  in 

Appendix  A. 

The  steady  state  solution  for  the  detailed  states  is 
presented  as  (Ref  11;253): 


P(xi,X2 xi.|)  = Cf  X ,)  f 2(  x^)  . . .f  Ki(  x^) 


(3) 


where  C is  the  normalizing  constant  which  causes  the  steady 
state  probabilities  to  sum  to  one  over  all  feasible  states. 
Each  function,  f , is  dependent  upon  the  type  of  queueing 
discipline  at  service  center  i.  For  a FCFS  service  center, 


j=1 


(4) 


For  a PS  service  center. 
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fi(xi)  = niifi  n 


r=l  1=1 


For  an  IS  service  center. 


^ i r i r 1 


irl 


Miri  1 


"’ikl- 


(5) 


R u . 

1 r 
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For  a Preemptive-Resume  LCFS  service  center. 


n ^ 

f . ( X . ) = 1 1 e . A . 
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The  specification  of  all  feasible  states  using  the 


detailed  state  description  is  extremely  time  consuming  since 


all  combinations  of  customers  must  be  considered  as  well  as 


the  permutations  of  customers  within  a particular  service 


center.  In  order  to  alleviate  this  problem,  Baskett,  et.  al. 


(Ref  11)  have  developed  an  aggregate  state  description 


The  aggregate  state  is  specified  by  S^=  ( y ■] , y 2 > • • • i Ym)  1 


where  y.=  ( n n n . „)  and  n.  is  the  number  of  class  r 

1 i1i2  iR  ir 


customers  at  service  center  i.  The  workload  parameter  for 


the  multiple  class  model  is  fJ  = ( n n n „)  . The  steady 

\ £1  A 


state  probabilities  are  described  by  the  product  form 


solution . 


P(S  = y ^ ,y^,  . . . ,yj^^)  = Cg  ^ ( y ^ ) g^C  yp)  . . . gj^,(  yj,^) 


where,  if  the  service  center  is  FCFS  then. 


II 

n1  p n . n . 

,&ir]  n/»i)  " 


if  the  service  center  is  PS  or  Preemptive-Resume  LCFS  then. 


,,, , = ,,.n  — 
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if  the  service  center  is  IS  then, 


h^h 


, , -.  n ^ M 


nj^!L“irJ 


(11) 


The  (Y^  ) are  determined  by  summing  the  T over  all  states 

which  have  n customers  at  center  i.  C in  Eq  (8)  is  the 
i 

normalization  constant  which  insures  that  the  sum  of  the 
probabilities  over  all  states  will  equal  one.  The 
normalization  constant  is  dependent  upon  the  nature  of  the 


workload,  therefore: 


C" ^ ( N ) = 


(y,  )tJ^(y2)-..G.,(y.,) 


(12) 


all  feasible 
states 


One  form  of  queue  dependent  service  distribution 
will  be  of  special  interest.  Let  be  the  rate  at 

which  one  server  at  the  i^  center  is  serving  class  r 
customers  when  3^  = ( n^  ^ , n^2  > • • • * r ) • ‘'hen  11^^^  ( ) = f ^ ( n^^ ) , 
where  fi(n-[)  is  an  arbitrary  positive  function  of  the  total 
number  of  customers  at  the  ith  center,  then  the  g^(S^)  in 
Eq  (8)  are  replaced  by: 


;*(S.)  = g.  (S.  )/n  f.  ( j) 


(13) 
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This  form  may  be  used  with  service  center  types  1, 
2,  or  4 (Ref  17;31)«  Other  forms  of  the  queue  dependent 
service  rates  are  presented  by  Daskett,  et.  al.  (Ref  11). 

III. 3 Consequences  o f Multiple  Classes 

An  important  assumption  of  the  solution  for  the 
single  class  model,  Appendix  A,  is  that  the  traffic 
intensity  for  a given  workload  remained  constant  in  the 
steady  state  condition.  This  assumption  permitted  the 
introduction  of  a new  variable,  Q ( n ^ n 2 , . • . n in  Eq  (32). 
The  traffic  intensity  is  defined  as  The  mean 
service  time  for  a server  at  the  i_th  service  center  is  given 
by  l/pj^,  and  the  mean  arrival  rate  of  customers  to  that 
center  is  For  a computer  system,  -jir  is  the  mean  time  the 
server  (resource)  at  center  i takes  to  satisfy  the 
requirements  of  a customer  if  the  full  capacity  of  the 
server  is  devoted  to  that  customer.  It  is  assumed  that  this 
time  does  not  change  with  the  workload  or  the  number  of 
customers  in  the  queue. 

when  more  than  one  class  is  considered,  the  traffic 
intensity  for  each  class  becomes,  X.  = a.  /p.  . The  mean 
service  time  and  mean  arrival  rate  are  now  specified  for 
each  class  at  the  iUi  service  center.  The  service  time  is 
still  dependent  upon  the  requirements  of  the  customer  and 
the  capacity  of  the  server.  However,  the  arrival  rate  may 
now  show  a dependence  upon  the  number  of  customers  from 
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i. 


other  dosses  present  at  the  queue  and  on  the  queueing 


discipline.  As  noted  earlier,  in  the  steady  state 
condition,  the  number  of  customers  of  a particular  class 
arriving  at  the  ith  service  center  must  equal  the  number  in 
that  class  departing  the  center.  When  processor  sharing  is 
used,  the  number  of  class  r departures  is  given  by: 

n. 

P(S^ 


where  p / is  the  fraction  of  the  server  capacity  devoted 
to  class  r customers.  An  identical  expression  may  be  used 
for  a fCFS  service  center.  In  this  case,  represents 
tlie  probability  that  a class  r customer  is  the  single 
customer  receiving  service  at  the  ith  center.  For  the 
general  case,  the  probability  that  class  r customers  are 
receiving  service  is  given  by: 


'Hr/ 


1 r 


'Hs/ 


1 s 


which  specifies  the  fraction  of  total  service  requirement 
devoted  to  class  r customers.  For  the  single  class  of 
customer  model  or  when  customers  have  identical  service  time 
requirements,  this  fraction  reduces  to  the  earlier  foi*m  used 
for  the  PS  service  center.  The  solutions  to  the  balance 
equations  are  valid  only  if  the  customers  at  a FCFS  service 
center  have  identical  service  time  reiiuirements . When  tlie 
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socvioe  time  reqvi  i reinents  at  a FCFS  center  are  not 
identical,  the  local  balance  equations  are  inconsistent 
(Ref  1 1 ;253) . 

Although  a closed  form  solution  is  not  possible 
using  the  technique  of  local  balance,  some  form  of 
approximation  may  be  possible.  Iterative  procedures  have 
been  used  to  give  approximate  results  for  a fixed  workload. 
Baskett,  et.  al.  (Ref  11)  used  the  predicted  arrival  rates 
to  compute  the  mean  service  time  and  transition 
probabilities  for  an  "average"  customer.  Chen  (Ref  18) 
used  an  iterative  approach  to  predict  performance  for  a 
model  with  queue  dependent  transition  probabilities. 

When  customers  at  a FCFS  sei'vice  center  do  not  have 
identical  service  time  distributions,  then  some  correction 
to  the  mean  service  time  of  each  class  is  required.  Consider 
two  classes  of  customers;  one  with  a large  mean  service 
time,  the  other  witli  a small  mean  service  time.  The  time  in 
the  queue  for  the  large  job  will  be  less  than  it  would  be  if 
all  the  service  time.",  were  large.  This  arises  from  the 
smaller  wait  time  required  while  a small  job  is  serviced. 
The  opposite  is  true  for  small  jobs.  They  will  tend  to  wait 
longer  in  the  queue  than  would  be  required  if  all  tlie  jobs 
wore  small.  Therefore  the  mean  sei-vioe  time  for  large  jobs 
needs  to  be  decreased  and  the  mean  service  time  for  small 
jobs  needs  to  be  increased  if  accurate  predictions  are  to  bo 
made.  The  amount  of  change  will  depend  ujion:  the  amount  of 
server  capacity  devoted  to  each  cliss  and  the  magnitude  of 
the  differences  in  service  requirements. 
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1 1 1 . 4 Service  T i t;i e A d strn ^ti^t 

The  applicability  oT  the  product  form  solution 
depends  in  large  part  upon  how  well  the  restrictions  on  the 
queueing  types  are  adhered  to.  In  particular,  the 
restriction  requiring  equal  service  time  distributions  for 
each  class  of  customer  at  a FCFS  service  center  will 
generally  not  be  met  at  a service  center  with  a Round-Robin 
scheduling  discipline.  Even  tliough  customers  may  be 
allocated  equal  time  quantum  for  processing  at  the  resource, 
the  amount  of  service  required  on  the  last  visit  to  the 
server  will  usually  vary  with  different  customers,  or  even 
with  customers  in  the  same  class.  Consider  a workload 
consisting  of  two  classes  of  customers;  tlie  first  class  of 
customers  have  a mean  service  requirement  of  7 units  while 
the  second  class  lias  a requirement  for  70  units  of  service. 
If  the  scl'ieduling  algorithm  assigns  a quantum  of  7 units  to 
eacli  class  of  customer  during  each  visit,  then  the  service 
distribution  requirements  for  a FCFS  service  center  are 
satisfied.  In  this  instance,  the  probability  tl'.at  a class  r 
customer  is  receiving  service  at  the  ith  service  center  is 
^’ir^'^i  ■ sclieduling  algorithm  assigns  a time 
quantum  of  70  units  to  both  classes,  tl'.en  this  simple 
expression  for  tlio  probability  of  receiving  service  is  no 
longer  applicable.  In  general,  the  probability  of  a class  r 
customer  being  the  one  receiving  service  at  a FCFS  single 
server  center  is  expressed  as  the  ratio  of  total  class  r 
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service  reciiiired  at  that  resource  divided  by  the  total 
service  required  by  all  classes  at  the  resource. 

The  involved  scheduling  procedures  implemented  on 
the  DECsystem- 1 0 introduce  additional  complications  which 
hinder  the  direct  application  of  a closed  queueing  network 
model.  This  is  especially  true  when  the  system  is  in  the 
dual  processor  configuration.  The  customers  in  PQ 1 , with 
their  small  time  quantum  allocation,  receive  preferential 
assignment  to  CPUO,  the  primary  processor.  Customers  in  PQ2 
are  scanned  first  for  allocation  of  CPUl,  the  secondary 
processor.  In  addition,  the  count  for  in-core-protect  time 
may  cause  a customer  to  be  requeued  even  tnougli  its 
allocated  processing  quantum  has  not  expired.  For  the 
synthetic  jobs  used  for  this  validation  study,  the  in-core- 
protect  time  expiration  caused  requeueing  only  during  the 
first  time  quantum  in  PQ2. 

The  service  requirements  of  a cusLotiier  in  P02  can  be 
viewed  as  consisting  of  the  portion  which  uses  the  full 
quantum  allocations  and  the  remainder  portion  for  the  last 
pass  through  the  round  robin  cycle.  Each  part  of  the  total 
Ps2  service  requirement  can  be  considered  as  an  individual 
customer  to  the  FCFS  service  center.  The  interaction  of  the 
partial  service  requirements  for  each  class  of  customer 
causes  some  portion  to  remain  in  the  queue  longer  than 
expected  if  true  Processor  Sharing  where  used.  Similarly, 
other  portions  may  complete  service  faster  than  expected  for 
a PS  service  center. 
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a small  remainder  portion  while  there  are  many  customers 
which  use  full  quantum  allocations.  This  point  is 
illust!'ated  by  a workload  consisting,  of  a sinp.le  cl.vss  of 
customers  witti  mean  service  requirement  of  1 units.  If  the 
service  time  quantum  is  set  at  JO  units,  then  tlie  customers 
will  have  to  process  throut^h  tlie  queue  a second  time  to 
receive  the  final  unit  of  service.  The  total  time  at  the 
service  center  may  be  much  lonj',er  tlvan  it  would  have  been  if 
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units.  The 
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present.ed  in 

Chapter  V. 

lir.h  Summary 

liaskett,  ot.  al.  (kef  M),  vievelopeti  a Si)lui.ion 
for  the  v'l»^sed  Oueueinp,  Netwi.>rk  Model  with  multiple  cl.isses 
of  customers.  Althoup,!)  i.ltis  solut.ion  permits  more 
flexibility  for  represent  inp,  an  actual  comput  ei'  sy;ilem, 
there  are  assumptions  required  t'or  the  soluiiin  which  a?'e 
not  totally  s.itisfied  in  real  world  sy;<tem;'.  Tliis  v'liapter 
has  presented  a discussion  of  the  r-olulion  vievcloped  by 
Haskett,  et.  al.  and  has  iiuiieated  the  physical  inpl  i c at.  ions 
of  the  assum['t  i on  s . In  addit.ion,  an  ex  i.en;’>  i ve  vi  i Si'iKO'' ton  of 
the  restrictions  on  i'CKM  service  ciuiters  is  presented.  In 
this  final  discussion,  the  limitations  of  t.he  ;?olulton,  as 
they  apply  to  the  liK.Cs ys tern- 1 0 , have  be<'n  emphas i /('d  . 
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Chapter  IV 

Computer  performance  Measu_r_es 

The  solution  of  the  equilibrium  state  probabilities 
presented  in  Chapter  III  permits  the  calculation  of  system 
performance  measures  such  as  marginal  queue  length 
distributions  for  each  queue,  resource  ut  i 1 i 7.at  ion  , and  mean 
response  time.  This  chapter  presents  computationally 
efficient  algorithms  for  computing  the  mean  arrival  rate, 

. , and  the  resource  utilization,  U : These  algorithms  are 
not  applicable  to  all  the  service  centers  discussed  in 
Chapter  III,  but  they  may  be  used  with  many  models  of 
interest  so  that  computation  requirements  may  bo  reduced. 

IV. 1 Normalization  Constant  Calculation 

In  order  to  obtain  the  equilibrium  state 
probabilities,  the  g^(S^)  must  be  calculated  for  all  the 
service  centers  and  all  states  of  the  center.  In  addition, 
the  normalization  constant  must  be  dolerminel  by  summing  the 
product  form  solution  over  all  feasible  states.  If  C(fl)  is 
computed  in  a straight-forward  manner,  t tie  amount  of 
computation  increases  proportional  to  tlio  total  number  of 
states  in  tlie  model.  If  the  network  is  closed  iiul  the 
customers  do  not  change  class,  and  '’very  custotner  visits 
every  service  center,  then  the  number  of  states  is  given  by 
(Ref  17; 43): 


4n 


( 


II 


"TTHT 


TV  ■Ti’’ 


For  models  of  interest,  the  straight-forward  computation  of 
C(N)  may  become  impractical.  Wong  (Ref  1?)  and  McKenzie 
(Ref  1)  present  a computationally  more  efficient  method 
for  determining  C(N).  Two  restrictions  are  placed  on  the 
use  of  this  method  (Ref  17;J|7): 

1)  Customers  do  not  change  class 

2)  Fvery  class  of  customer  visits  every 
service  center 

Adjustments  may  be  made  in  ttie  network's  state  description 
to  overcome  the  first  restriction,  while  relatively  minor 
clianges  in  the  comput.ational  method  will  permit  customer 
classes  which  do  not  visit  every  service  center.  The  program 
developed  by  Mckenzie  (Ref  1)  requires  tliat  botli 
restrictions  be  met.  For  tlie  purpose  of  this  thesis,  tliese 
restrictions  were  not  unreasonable.  More  detailed  models, 
those  modeling  a greater  variety  of  system  resources,  may 
require  that  class  changjes  be  permitted  or  that  certain 
classes  do  not  visit  all  service  centers.  These  models  will 
require  that  the  program  be  modified  to  include  the 
alternate  procedures  suggesteti  by  Wong  (Ref  17;b0). 

[V.2  Qnpue  I.engtli  Calculation 

After  the  equilibrium  state  probabilities  and 

normalization  constant  have  been  calculated,  t(,i  .’’argina' 
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queue  length  distributions  and  expected  number  of  customers 
can  be  computed  for  each  resource  and  each  class  of 
customer.  The  probability  that  there  are  n customers  at  the 
ith  service  center  is  given  by  (Ref  17;54): 


P.(n:)  = 2 Pi(S.) 
all  S 
such  that 

n +n  +...+n  =n 
i1  i2  iR  i 


(15) 


The  probability  that  there  are  n^^  class  r customers  at  the 
ith  service  center  is  given  by: 
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(16) 


The  expected  number  of  total  customers,  n^,  and  expected 

number  of  class  r customers,  n.  , at  tiie  ith  service  center 

ir’ 

are  easily  computed  from  Eqs  (15)  and  (16). 


I V . 3 Re  sour  c^e  U t i l_i  z a ^ijJ  n C a l_c  u 1 a 1 1 

For  a single  server  service  center,  resource 

utilization,  IJ  . , is  the  fraction  of  a unit  time  interval 
ir  ’ 

that  the  server  at  ith  center  is  busy  serving  class  r 
customers.  When  the  service  cetiter  is  type  1,  2,  or  '1, 
is  computed  by: 
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where  =0  when  and  n^  are  both  zero.  This  ratio  is 
the  fraction  of  server  capacity  received  by  class  r 
customers  when  the  iUi  center  has  a PS  discipline.  For  FCFS 
or  Preemptive-Resume  LCFS  centers,  the  ratio  is  the 
probability  that  a class  r customer  is  getting  full  capacity 
of  the  i_th  service  center.  This  assumes  that  the  service 
time  distributions  meet  the  restrictions  described  in 
Chapter  III. 

If  the  service  center  is  single  server  and  has  a 
constant  service  rate,  then  a simplier  expression  for 
resource  utilization  is  possible.  Resource  utilization  may 
be  expressed  as  (Ref  17;55): 


ftong  (Ref  17)  expands  on  this  expression  after  he  proves 
that  the  arrival  rate  is  given  by: 


C(N) 


(19) 


where  N *”=  ( N ^ iJ  2i  • • . N 1 , . . , N j^)  , which  permits  the 
utilization  to  be  expressed  by: 


When  service  center  i contains  a.  identical 

i 

resources,  the  utilization,  is  the  average  measure  of 

tlie  utilization  of  the  individual  resources.  Therefore: 

Uif,  = (21) 

^i^r 

If  the  service  center  has  a PS  discipline,  then  the 
following  expression  holds,  even  for  multiple  classes  of 
customers : 


C(N) 

C(N-'') 


(22) 


For  service  centers  of  type  1 or  4,  this  expression  does  not 
apply.  The  utilization,  must  be  computed  by  Fq  (17). 
The  straight-forward  approach  must  also  be  used  when  the 
servers  have  state  dependent  service  rates.  The  processors 
on  the  DECsystem-10  at  the  Air  Force  Avionics  Laboratory  are 
not  identical.  The  operating  system  overhead  on  the  primary 
processor  is  consistently  higher  than  that  of  the  secondary 
processor.  In  addition,  the  secondary  must  stop  processing 
a customer  which  initiates  an  I/O  request  since  only  tlie 
primary  has  the  normal  I/O  capability.  Because  of  these 
features,  the  model  for  the  pr imary /secondary  con f igur ation 
should  be  state  dependent  in  order  to  accurately  reflect  the 


physical  system. 
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The  user  terminals  in  a time-sharing  computer  system 

are  represented  by  an  Infinite  Server  (IS)  service  center. 

This  type  of  service  center  has  an  expression  for  resource 

utilization  similar  to  the  expressions  for  the  other  service 

center  types.  In  an  IS  service  center,  there  are  at  least 

N.,  servers  dedicated  to  the  N class  r customers  in  the 
f r 

system.  If  the  IS  service  center  were  represented  as  having 

identical  servers  then  the  utilization  would  be 

A.  /N  M.  . However  it  is  notationally  more  convenient  to 
1 r r 1 r 

express  the  utilisation  as  (Ref  17;58): 
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where  will  be  the  mean  number  of  class  r customers  at 
the  i^  center  when  it  is  an  IS  service  center.  Then  the 
utilization  can  be  expressed  as: 


These  expression  for  a and  U allow  for  more 

i r ir 

efficient  computation  of  the  moan  arrival  rate  and  resource 
utilization.  The  values  of  C(N“^)  are  derived  in  t’r.e 
iterative  process  used  by  McKenzie  (Ref  1)  to  calculate 
the  normalization  constant.  Wicn  tiie  service  center  is  of 
the  type  discussed  above,  the  values  of  C(N~*')  can  be 
retained  and  the  less  efficient  computations  required  by 
Eq  (17)  can  be  avoided. 
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IV.  4 Response  Tiine  Cjilculation 

Ivhen  has  been  calculated,  the  response  time  for 

class  r customers  can  be  determined.  Response  time  is  the 

amount  of  time  a customer  spends  in  the  network  outside  of 

the  IS  service  center  representing  the  terminals.  This  is 

shown  in  Figure  6.  The  mean  number  of  class  r customers 

"thinking"  at  the  terminals  is  given  by  The  number  of 

class  r customers  in  the  remainder  of  the  network  is  N -n.  . 

r 1 r 

From  Little's  formula  (Ref  19),  the  mean  response  time  is 
given  by : 


ir 


N p-n 

^r 


(25) 
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IV.  5 

The  number  of  possible  states  for  a Closed  Queueing 
Network  Model  increases  considerably  with  the  addition  of 
service  centers  or  the  spec i Cic ation  of  multiple  classes. 
The  direct  computation  of  the  normal  i zati(5n  constant  of  the 
product  form  solution  and  resource  utilization  requires  that 
specific  values  be  summed  over  all  possible  states.  In  order 
to  reduce  the  computational  requirements  for  solution  of 
these  models,  alternate  methods  have  been  presented  to 
determine  the  normalization  constant,  mean  arrival  rate  and 
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resource  utilization.  These  methods  do  require  that 
additional  assumptions  be  made  concerning  the  model. 
However,  they  are  applicable  to  most  models  of  interest. 
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£hapte£  V 
Model  Validation 


The  concept  of  model  validation  is  not  easily 
defined.  However,  "the  validation  ought  to  establish  some 
degree  of  confidence  that  the  model  produces  results  which 
sufficiently  mirror  the  significant  attributes  of  the  system 
being  modeled  (Ref  21; 139)."  The  time  limitations  of  this 
work  prohibit  the  use  of  detailed  hypothesis  testing 
techniques  which  should  be  used  to  make  a final  model 
selection.  Since  this  report  is  primarily  concerned  with 
performance  for  time-sharing  users,  value  judgments  based  on 
response  time  errors  will  be  the  basis  for  evaluating  the 
quality  of  a model. 

It  is  not  possible  to  validate  all  combinations  of 
system  characteristics  with  absolute  certainty.  In  any  case, 
many  of  the  resource  configuration  and  workload  combinations 
would  be  of  little  value  to  understanding  tiie  real  system. 
The  validation  effort  should  be  directed  towards  common 
system  con f igurat ions  and  reasonable  workload 
characterizations.  The  model  validation  effort  investigates 
the  application  of  several  Closed  Queueing  Network  Models  to 
the  DECsystom-10  dual  processor  system.  In  addition,  the 
single  processor  configuration  is  examined  as  well  as  the 
application  of  a model  to  represent  the  system  when  swapping 
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occurs . 


’;  The  validation  of  the  Closed  Queueing  Network  Model 

is  performed  with  the  use  of  a synthetic  workload.  The 
workload  is  developed  using  the  SCRIPT  program  to  run 
; duplicate  copies  of  selected  synthetic  jobs.  Each  job  cycles 

through  a predetermined  sequence  of  subroutine  calls.  The 
cycle  starts  with  the  job  in  a sleep  state  which  is  used  to 
model  teletype  think  time.  After  the  sleep  state,  a CPU 
subroutine  is  called  a specified  number  of  times.  Finally, 
the  job  makes  a specified  number  of  calls  to  a disk  I/O 
i subroutine  which  writes  a single  block  of  data  onto  a 

designated  RP03  disk  file.  'when  tlie  disk  I/O  is  completed, 
the  job  returns  to  the  sleep  state  which  represents  user 
; think  time.  Numerous  combinations  of  processing  and  I/O 

; requirements  are  possible  with  synthetic  jobs  of  this  form. 

[ One  restriction  of  the  synthetic  jobs  which  is  not  required 

i of  real  jobs  is  that  processing  and  disk  I/O  are  sequential, 

i Although  processing  and  I/O  may  overlap  when  different  jobs 

■'  are  involved,  processing  and  I/O  within  a single  job  do  not 

' overlap.  This  restriction  on  the  synthetic  workload 

; coincides  with  the  queueing  model  restriction  prohibiting  a 

* 

customer  from  using  two  resources  at  one  time. 

The  overhead  introduced  by  the  Operating  System  is 
represented  by  increasing  the  processing  service 
■ requirements  by  the  amount  of  service  time  which  would  be 

I 

lost  to  overhead.  Because  the  Operating  System  overhead  is 
not  equally  divided  between  the  two  processors,  they  will 
f not  be  considered  as  identical  servers.  The  use  of  state 

j 
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dependent  service  rates  does  permit  realistic  modeling  of 


*. 


this  aspect  of  the  processors.  Customers  in  PQ2  are 
considered  first  for  assignment  to  the  secondary  processor. 
Therefore  the  overhead  of  CPU1,  approximately  5%,  is  used  to 
determine  the  increase  in  service  time  requirements  of  each 
customer  in  PQ2.  The  overhead  of  the  primary  processor  is 
much  higher,  approximately  12?.  Because  of  the  larger 
overhead  on  CPUO,  the  service  rate  is  multiplied  by  1.93 
when  two  or  more  customers  are  present  in  the  queue.  The 
service  time  requirement  of  customers  in  PQ1  is  increased  by 
an  amount  determined  by  CPUO  overhead. 

This  validation  study  considers  two  basic  questions 
raised  by  Giammo  (Ref  21): 


1)  Do  the  service  time  assumptions 
introduce  a basic  distortion? 

2)  Can  the  logic  of  job  processing  be 
sufficiently  modeled  as  a network  of  queues? 


The  validation  effort  considers  numerous  workload 
combinations  and  the  basic  processor  configurations 
available  on  the  DECsystem- 1 0 . Although  disk  I/O  is  included 
in  the  workload  requirements,  the  main  emphasis  of  this  work 
is  directed  towards  accurately  modeling  the  processing 
capabilities  of  the  system.  The  involved  scheduling  and 
queueing  procedures  associated  with  the  processors  justify 
detailed  analysis  of  this  area.  A total  of  21  experiments 
were  conducted  to  validate  tne  models  examined  in  this 
report . 
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Knch  validjtioii  oxpefimoiit  coiisidofu  a specific 
model  conf  igiifation  and  combinations  of  two  or  three 
synthetic  workloads.  The  workloads  consist  of  various 
combinations  of  customers  from  one  or  two  classes.  Kach 
customer  class  has  predetermined  processini;  and  disk  I/O 
r equ  i I'cments  . The  individual  experiments  will  analyze  the 
effects  of  different  numbers  of  customers  from  each  class 
and  ch.-ini'es  in  tlie  processini',  requirements  of  the  classes. 
For  the  jobs  executed  in  this  validation  work,  the 
nFCsystem-lO  run  queue  is  considered  in  two  portions,  PQ1 
and  PQP.  oince  Jobs  in  eac^h  iiucue  are  handled  differently  by 
the  schedulei',  a natural  viuestion  concerns  how  processor 
viueueini;  should  be  represented.  The  models  evaluated  in  tliis 
repoi't  use  two  seperate  service  centers  in  one  ca:3e  and  in 
the  otlier,  a sitq^le  service  center  to  represent  the  queues 
wliich  foraii  at  the  processor  resource. 

V.1  Two  Node  Proct^ssor  Model 

Fxperiments  1 throuph  6 i nv  est  i p,  at  ed  the  effect  of 
ilifforent  workloads  on  a specific  moilel  of  the  IMiCsystem- 1 0 
in  the  dual  processor  con  f ip.urat  ion  . Facli  experiment 
considers  combinations  of  the  synthetic  jobs  with 
(iredi>term  i netl  process  irq;  ntid  I /t>  requ  i rem<'n  ts  . The  model 
con  f i j^ur  at  ion  u..ea  for  tliis  series  of  experiments  is  l.he 
multiple  r.erver-  model  used  by  Mckenzie  (lief  ll.  This 
model,  shown  in  Fip.ure  Y,  represents  tlio  telei.ype  by  an 
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igure  7.  Two  ;iode  .Processor  Xodel 


Infinite  Server  service  center.  The  processors  are 
represented  by  a serial  combination  of  two  service  centers. 
A single  server  FCFS  center  services  customers  in  PQ1  as 
they  depart  the  teletype  node.  Customers  in  PQ2  are 
serviced  by  a multiple  server  Processor  Sharing  center.  Disk 
I/O  is  re presented  by  a single  server  FCFS  center  with 
feedback  to  represent  multiple  disk  calls.  The  service 
times  and  transition  probabilities  in  Figure  7 are  those 
used  for  experiment  6. 

V. 1 . 1 Experiment  2 

Experiment  1 investigates  the  application  of  the 
model  shown  in  Figure  7 for  a workload  which  consists  of  a 
single  class  of  job  with  the  following  characteristics: 

Think  Time  = 300.0  jiffies 
Processor  Service  = 7.6  jiffies 
Disk  I/O  Calls  = 3 


The  results  of  this  experiment  are  shown  in  Table  I,  where 
the  workload  consisted  of  20  jobs  for  part  1 and  10  jobs  for 
part  2.  The  workload  closely  matches  the  service  time 
distribution  requirement  placed  on  a FCFS  service  center. 
Although  the  model  uses  a Processor  Sharing  service  center 
to  represent  the  processing  in  PQ2,  tiie  real  system's  round- 
robin  scheduling  .more  closely  resembles  a FCFS  center  with 
feedback,  as  discussed  in  Chapter  III.  This  concept  of  the 
actual  system  will  permit  a better  understanding  of  the 
results  of  this  validation  study.  Although  the  equal  service 
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times  of  the  jobs  would  seem  to  exactly  satisfy  the 
requirements  for  a FCFS  service  center,  the  sequencing  of 
jobs  through  PQ1  and  PQ2  may  cause  the  jobs  to  have 
differing  service  requirements  depending  on  which  cycle  of 
the  Round-Robin  sequence  they  are  in.  In  this  example  each 
job  will  receive  3 jiffies  of  service  while  in  PQ1  and  the 
remaining  4.6  jiffies  while  in  PQ2;  service  time 
compensation  for  overhead  will  be  ignored  for  this 
discussion.  The  minimum  PQ2  time  quantum  is  larger  than  the 
service  requirement  remaining  after  the  job  leaves  Pyi. 
Therefore,  each  job  will  receive  two  relatively  equal  bursts 
of  service  from  the  processors,  one  of  3 jiffies  and  the 
other  of  4.6  jiffies.  This  factor  contributes  to  the 
extremely  low  response  time  error  shown  for  both  workloads. 

V.  1 . 2 Kxperiiiient  2 

The  workload  for  experiment  2 consists  of  a single 
class  of  customers  with  the  following  cliaracter  istics : 

Think  time  = 300.0  jiffies 
Processor  Service  = 30.6  jiffies 
Disk  I/O  Calls  = 10 

The  results  of  this  experiment  are  shown  in  Table  II  for  20 
jobs  and  10  jobs  in  part  1 and  part  2 respectively.  The 
response  time  error  has  increased  for  both  workload  levels. 
The  processing  requirement  of  the  jobs  in  this  experiment 
produce  processing  cycles  of  differing  length,  even  though 
the  total  service  times  of  the  jobs  are  identical.  The  PQ1 
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Table  II 
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time  slice  is  3 jiffies,  while  the  remainder  received  in  PQ2 


is  27.6  jiffies.  The  assumption  of  the  model  is  that  all 
jobs  either  receive  equal  service  of  the  resource,  for 
Processor  Sharing,  or  that  the  probability  of  receiving 
service  is  equal  for  all  jobs,  for  a FCFS  resource.  The 
differing  quantum  allocations  prevent  this  assumption  from 
being  met.  However,  the  priority  assignment  of  jobs  in  PQ1 
to  CPUO,  and  jobs  in  PQ2  to  CPU1  does  compensate  somewhat  by 
providing  for  jobs  in  PQ1  a more  equal  probability  of 
receiving  service.  McKenzie  (Ref  1)  discusses  one 
procedure  to  help  alleviate  this  problem.  In  effect,  the 
time  that  jobs  in  PQ2  wait  for  jobs  in  Pyi  is  added  to  the 
total  service  time.  However,  this  procedure  does  not  address 
all  aspects  of  the  problem.  The  service  time  needs  to  be 
decreased  by  some  amount  because  jobs  in  PQ2  on  the  real 
world  system  must  wait  for  a job  with  a time  slice  of  only 
27.6  rather  than  the  30.6  assumed  by  the  model.  Offsetting 
this  adjustment  is  the  added  time  that  jobs  in  PQ1,  with  a 3 
jiffy  time  slice,  must  wait  for  jobs  from  PQ2  which  may  liave 
been  assigned  to  both  processors.  The  effect  of  each  of 
these  factors  is  difficult  to  determine,  but  it  would  seem 
to  depend  upon  how  busy  the  processors  are  u'ith  each  type  of 
job.  For  this  particular  experiment,  the  model  predictions 
are  less  than  measured  when  the  processors  arc  very  busy; 
while  the  predictions  are  higher  than  measured  i.'hen  there  is 
low  processor  utilization. 

A final  consideration  which  applies  to  all 
experiments  which  display  relatively  low  processor 
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utilization  is  that  the  deterministic  nature  of  the 

* synthetic  jobs  does  not  tjuarantee  that  a true  steady  state 

condition  is  obtained.  In  fact,  observations  of  the 

processor  control  display  during  these  experiments  indicated 
that  jobs  tended  to  process  in  "bunches'*  when  the  processor 
utilization  was  approximately  85%  or  lower.  This 


observation  was  further  substantiated  by  the  METER  reports 
which  indicated  that  the  response  time  of  the  synthetic  jobs 
tended  to  be  in  one  of  two  clusters  under  this  workload 
condition.  The  measured  response  time  for  each  cycle  of  the 
synthetic  job  in  experiment  2.2  fell  in  a groups  centered  at 
approximately  135  jiffies  and  75  jiffies.  If  the  conditions 
of  the  model  were  completely  satisfied,  the  measured 
response  times  would  be  more  closely  distributed  about  the 
overall  mean. 

V.1.3  Experiment  3 

The  workload  for  experiment  3 also  consists  of  a 
single  class  of  customer  with  the  following  characteristics: 

Think  Time  = 300.0  jiffies 

Processor  Service  = 83.0  jiffies 

Disk  I/O  Calls  = 3 

The  results  of  this  experiment  for  20  and  10  jobs  are  shown 
in  Table  III  part  1 and  part  2 respectively.  The  response 


time  errors  are  similar  to  those  shown  for  experiment  2.  For 
this  experiment,  the  longer  service  time  required  at  the 
processor  permits  each  job  to  receive  more  time  quanta  while 
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in  PQ2.  The  lont'er  service  time  also  ctiuses  more  jobs  to 
wait  in  PQ2  and  increases  the  probability  that  jobs  from  PQ2 
will  be  using  both  processors.  Because  of  this,  jobs  coming 
out  of  the  sleep  state  will  more  likely  be  required  to  wait 
while  a job  from  PQ2  completes  service  on  CPUO,  but  the 
overall  processor  service  will  be  closer  to  processor 
sharing . 

V.1.4  Expej'iment  4 

The  job  mix  for  experiment  4 consists  of  two  classes 
of  jobs  with  the  following  characteristics: 


Think  Time  = 

Processor  Service 

Class  1 = 
Class  2 = 
Disk  I/O  Calls 

Class  1 = 
Class  2 = 


300.0  jiffies 

38.0  jiffies 

76.0  jiffies 

10 

3 


The  results  of  this  experiment  are  shown  in  Table  IV  parts 
1,  2 and  3 for  job  mixes  of  10  each,  8 each,  and  5 each 
respectively.  The  use  of  this  more  realistic  workload 
provides  a clearer  picture  of  the  value  of  the  Closed 
Queueing  Network  Model.  kitli  this  workload,  the  actual 
computer  system  has  a variety  of  customers  in  differing 
stages  of  completion  in  tlie  Pound-Robin  processing  cycle. 
This  provides  an  even  greater  departure  from  the  equal 
service  rate  requirement  for  a true  FCf'S  service  center.  The 
good  error  rate  for  the  mix  of  5 jobs  of  each  type  is 
believed  to  be  the  result  of  the  failure  to  reach  a true 
steady  state  condition  rather  than  improved  model  accuracy. 
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V . 1 . 5 Experiment  5 

Two  classes  of  jobs  were  also  used  for  experiment  5. 


These  jobs  have  the  following  characteristics: 


Think  Time 
Processor  Service 
Class  1 
Class  2 

Disk  I/O  Calls 

Class  1 
Class  2 


300.0  jiffies 

30.6  jiffies 

83. 6 jiffies 

10 

3 


The  results  of  this  experiment  for  job  mixes  of  10  each,  8 
each,  and  5 each  are  shown  in  Table  V parts  1,  2 and  3.  A 
slight  improvement  in  response  time  error  is  indicated  over 
similar  job  mixes  in  experiment  4.  The  time  quantum 
allocated  for  both  classes  of  jobs  is  more  fully  utilized 
than  was  true  for  experiment  4.  For  example,  type  2 
customers  in  experiment  4 had  a PQ2  remainder  cycle  of  5 
jiffies  which  competed  for  the  processors  with  other  jobs 
requiring  a full  30  jiffy  time  slice.  For  this  experiment, 
the  PQ2  service  requirement  of  type  1 jobs,  27.6  jiffies,  is 
satisfied  in  one  visit  to  the  processor.  The  changes  which 
produced  the  slight  improvement  for  response  time  error  in 
experiments  5.1  and  5.2,  over  similar  parts  of  experiment  4, 
had  a negative  effect  on  the  apparent  accuracy  of  the  model 
for  the  mix  of  5 jobs  of  each  class. 


V.1.6  Expej^iment  6 

The  final  experiment  using  the  model  represented  in 
Figure  7,  to  predict  the  performance  of  the  DF.Csystem- 1 0 in 
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Experiment  5 Results 
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the  dual  processor  configuration  also  uses  a job  mix  of  two 

class  types.  The  characteristics  of  the  jobs  used  in 

experiment  6 are  as  follows: 

Think  Time  = 300.0  jiffies 

Processor  Service 

Class  1 = 38.0  jiffies 

Class  2 = 7.6  jiffies 

Disk  I/O  Calls 

Class  1 = 10 

Class  2 = 3 

The  results  of  this  experiment  are  found  in  Table  VI  parts 
1,  2 and  3 for  job  mixes  of  10  each,  8 each,  and  5 each 
respectively.  This  experiment  considers  jobs  which  place 
light  demands  on  the  system.  Only  one  full  PQ2  time  quantum 
is  used,  that  being  for  type  1 jobs.  The  remaining  portions 
of  service  which  the  jobs  receive  in  the  Round-Robin 
processing  cycle  are  relatively  equal.  The  fact  that  most 
of  the  processing  cycles  are  nearly  equal  is  reflected  in 
the  good  response  time  errors. 

V . 2 0^  ilode  Processor  Model 

For  experiments  7 through  12,  the  DECsystem-10  dual 
processors  are  modeled  with  a single  service  center  using 
the  Processor  Sharing  queueing  discipline.  The  total  CPU 
service  time  requirement,  PQ1  time  plus  PQ2  time,  is  assumed 
to  be  received  at  this  service  center,  rather  than  in  stages 
as  in  experiments  1 through  6.  Figure  8 shows  this  model 
representation  of  the  DECsystem-10.  The  service  times  and 
transition  probabilities  in  the  figure  are  those  used  for 
experiment  12. 
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This  model  does  not  provide  for  the  priority  service 
which  jobs  receive  while  in  PQ1.  The  failure  to  provide 
priority  service  for  jobs  in  PQ1  tends  to  increase  predicted 
response  time  because,  in  the  model,  these  jobs  share  the 
total  capacity  of  the  processors  with  jobs  from  PQ2.  This 
same  feature  could  also  cause  some  decrease  in  predicted 
response  time  since  the  model  allows  jobs  in  PQ2  to  be 
processed  by  either  processor,  giving  equal  capacity  to  time 
in  PQI  and  PQ2.  This  feature  could  provide  jobs  at  the  CPU 
service  center  with  more  resource  service  than  would  be 
found  in  the  real  world  system.  Combining  the  two  run  queues 
will  have  different  effects  on  different  jobs  processed 
through  the  model.  The  direction  and  magnitude  of  error  will 
depend  upon  such  features  as  number  of  customers  receiving 
full  PQ2  time  quantum  and  the  number  of  customers  receiving 
service  while  in  PQ1.  The  results  of  experiments  7 through 
12  indicate  this  dependency  upon  the  type  of  workload. 
However,  these  results  also  show  that  this  simplied  system 
representation  is  approximately  as  accurate  as  the  model 
which  used  two  nodes  to  represent  the  processors. 


V.2.1  Z 

The  results  of  experiment  7 are  shown  in  Table  VII. 
For  this  experiment,  the  workload  characteristics  wore  the 
same  as  used  for  experiment  1.  Although  the  response  time 
errors  have  increased  for  both  workloads,  the  size  of  the 
error  is  still  relatively  small  when  compared  with  results 
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of  experiments  involving  more  complex  workloads.  The  fact 
that  the  predicted  response  times  decreased  further  below 
the  measured  values  indicates  that,  for  this  experiment,  the 
assumption  of  the  model  that  all  jobs  receive  an  equal  share 
of  capacity  from  the  processors  is  too  optimistic.  The 
measured  values  indicate  that  there  is  more  conflict  for 
processor  service  than  is  predicted  by  the  model. 

V . 2 . 2 Experiment  8 

The  results  for  experiment  8 are  shown  in  Table 
VIII.  This  experiment  used  the  same  workload  inputs  as 
experiment  2.  The  increase  in  predicted  response  times 
indicates  that  the  modeled  system  has  become  more  congested. 
This  results  from  the  increased  service  time  which  causes 
jobs  to  remain  in  the  run  queue  longer  and  therefore 
increases  the  number  of  jobs  in  the  queue  for  the 
processors.  Since  one  of  the  response  time  errors  decreased 
while  the  other  increased,  it  is  difficult  to  make  a 
comparison  between  the  two  models.  It  should  bo  noted  again 
however,  that  the  measured  data  for  light  workloads  did  not 
appear  to  be  from  a true  steady  state  condition.  Because  of 
this,  all  experiments  which  have  low  CPU  utilizations  are  in 
question.  These  experiments  may  still  be  of  interest  in  that 
they  establish  some  measure  of  accuracy  for  this  level  of 
workload  when  a steady  state  condition  is  not  achieved. 
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V .2.3  Experiment  9 

The  results  of  experiment  9 are  shown  in  Table  IX. 
The  workloads  used  in  this  experiment  were  the  same  as  used 
for  experiment  3.  The  response  time  errors  show  changes  very 
similar  to  those  of  experiment  8.  Since  the  total  CPU 
service  time  is  larger  than  used  for  experiment  8,  there  are 
more  jobs  in  the  run  queue.  In  the  model  this  results  in  a 
lower  percentage  of  server  capacity  for  jobs  waiting  at  the 
processor  service  center. 

V . 2 . 4 Experiment  1 0 

The  results  of  experiment  10  are  shown  in  Table  X. 
This  experiment  used  the  same  workloads  as  used  for 
experiment  4.  For  the  workloads  which  cause  heavy  CPU 
utilization,  the  response  time  errors  improved  for  both 
classes  of  customers.  For  this  experiment,  the  simplified 
model  appears  to  be  at  least  as  accurate  as  the  model  used 
in  experiment  4. 

V . 2 . 5 Experiment  1 1 

Table  XI  shows  the  results  for  experiment  11.  The 
workloads  for  this  experiment  were  the  same  as  used  for 
experiment  5.  The  response  time  errors  for  this  experiment 
show  the  same  trend  as  noted  for  experiment  10.  The  use  of 
one  Processor  Sharing  node  is  at  least  as  accurate  a 
representation  of  the  DECsystem-10  as  the  dual  node 
representation  used  for  experiment  5. 
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V . 2 . 6 Experiment  1 2 

The  results  of  experiment  12  are  shown  in  Table  XII. 
The  workloads  for  this  experiment  are  the  same  as  used  for 
experiment  6.  For  this  experiment,  the  one  processor  node 
representation  was  again  at  least  as  accurate  as  the  model 
with  two  processor  nodes. 

V . 3 Appl ication  of  Service  Time  Adjustment 

In  Chapter  III,  the  problem  of  unequal  service 
requirements  for  a Round-Robin  queueing  discipline  was 
discussed.  Eq  (14)  was  presented  as  a method  for  adjusting 
the  service  times  to  compensate  for  differences  between  the 
real  world  system  and  the  assumptions  imposed  by  the  model. 

The  application  of  Eq  (14)  requires  that  at  least 
one  model  run  be  made  without  the  service  time  adjustments. 
This  initial  model  run  permits  the  CPU  utilization  to  be 
predicted.  The  service  time  of  each  class  of  customer  is 
then  divided  into  the  portions  of  service  received  while  in 
PQ1  and  for  each  cycle  through  the  Round-Robin  scheduling  of 
PQ2.  Adjustments  are  made  to  each  portion  of  service 
received  in  PQ2  and  the  new  service  time  requirements 
determined  by  adding  the  individual  parts.  Additional 
service  also  is  determined  for  the  overhead  of  the  Operating 
System.  Once  computed,  the  adjusted  service  time 
requirements  can  be  applied  to  a model.  This  procedure  was 
applied  to  the  model  and  workloads  used  for  experiments  1 
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Experiment  12  Results 
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through  6.  The  service  times  and  transition  probabilities  in 
Figure  9 are  those  used  for  experiment  18.  The  results  of 
this  series  of  experiments  are  found  in  Table  XIII  through 
XVIII. 

For  most  applications  of  the  service  time 
adjustment,  the  response  time  errors  remained  about  the  same 
or  slightly  higher  than  those  of  the  corresponding 
experiments  1 through  6.  However,  the  application  of  service 
time  adjustments  to  the  workload  of  experiment  '4  caused  the 
response  time  errors  to  decrease  by  a significant  amount 
when  the  processors  were  heavily  utilized.  These  results 
appear  in  Table  XVI.  The  service  time  requirement  of  the 
jobs  in  that  experiment  were  such  that  that  the  last  cycle 
through  the  Round-Robin  PQ2  queue  required  a small  amount  of 
service  when  compared  with  the  full  quantum  allocated  to  the 
jobs  in  the  queue.  As  a result,  jobs  on  the  real  world 
system  had  to  wait  in  the  queue  to  receive  a small  amount  of 
service  from  the  processors.  The  use  of  the  service  time 
adjustment  is  an  attempt  to  correct  for  this  additional  wait 
not  incorporated  in  the  model.  The  procedure  shows  good 
improvement  in  response  time  errors  for  jobs  which  had  a 
small  service  time  requirement  for  the  final  pass  through 
the  Round-Robin  cycle  of  PQ2.  There  is  some  negative  impact 
on  other  job  mixes  but  the  effect  on  the  workloads  used  in 
this  report  were  relatively  small. 
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Experiment  14  Results 
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V.4  Single  Processor  Model 

The  service  time  adjustment  procedure  was  also 
applied  to  a model  representing  the  DKCsystem-10  with  only 
one  processor.  This  representation  is  shown  in  Figure  10. 
The  service  times  and  transition  probabilities  in  this 
figure  are  those  used  for  experiment  20. 

Time-sharing  users  may  have  only  one  CPU  available 
to  them  if  the  secondary  processor  malfunctions  or  if  the 
real  time  simulation  programs  are  being  run  at  high 
priority.  Experimetits  19  and  20  were  conducted  with  the 
secondary  CPU  switched  off-line.  The  model  was  evaluated 
once  for  the  predetermined  service  time  requirements,  plus 
overhead,  and  again  for  the  adjusted  service  time 
requirement,  calculated  with  Eq  ( 1 'I ) . Experiment  19 
considered  a sitigle  class  of  customer  while  experiment  20 
considered  two  customer  classes. 

V.4.1  Experiment  2*? 

The  workload  for  experiment  19  considered  a single 

class  of  10  jobs  with  the  following  characteristics: 

Think  Time  = 300.0  jiffies 
Processor  Service  = 83*b  jiffies 
Disk  I/O  Calls  = 3 

The  results  of  the  experiment  are  shown  in  Table  XIX.  Part  1 
presents  the  results  using  the  required  processor  service 
time  plus  the  overhead  correction.  This  direct  application 
of  the  model  to  the  single  processor  configuration  produces 
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results  romsrkably  close  to  measured  values.  This  is 
partially  accounted  for  by  the  relatively  equal  amounts  of 
service  required  each  time  throujih  the  Round-Robin  queue  for 
rQ2.  Each  customer  should  receive  two  allocations  of  30 
jiffies  and  have  a final  requirement  for  20.6  jiffies  while 
in  PvjP.  However,  since  each  job  receives  3 jiffies  of 
service  wlUle  in  POl,  ther-e  should  be  some  conflict  for  the 
use  of  the  single  processor  on  the  real  world  system.  In 
particular,  even  ttiough  jobs  in  Ps'l  have  priority  over  P02 
jobs,  the  jobs  enteritik;  POl  may  still  have  to  wait  until  the 
existing;  jobs  at  the  resource  completes  its  sei'vice  quantum. 
Since  the  processor  may  be  servicing  a job  from  PQ2,  the 
wait  time  could  be  as  long  as  30  jiffies.  The  jobs  in  t!ie 
real  world  system  could  then  be  expected  to  experience 
additional  waiting  not  included  in  the  model  solution.  The 
small  error  value  in  this  experiment  may  not  be  influenced 
so  much  by  the  accuracy  of  the  movlel  as  by  other  offsetting 
factors,  such  as  those  to  be  discussed  in  experiment  20. 

The  results  in  part  2 were  obtained  after  a 
correction  was  made  for  the  differences  in  service 
requirements  on  different  passes  throug,!!  t!u'  Round-Robin 
queue  of  P02.  There  was  some  improvement  for  nearly  all  the 
predicted  values.  Since  the  jobs  in  this  experiment  had  a 
large  remaining  service  requirement  on  tlie  final  pass 
through  the  Round-Robin  queue,  the  service  time  adj\istments 
were  small,  as  indicated  by  the  relatively  small  changes  in 
the  predicted  values.  The  changes  were  however,  in  a 
direction  which  produced  smeller  errors. 


l!UHlJll«mj!PH.iimii|_i.  ..i 


V . . 2 P'?'" 

The  results  for  experiment  20  are  shown  in  Table  XX. 
The  experiment  used  two  10  job  classes  with  the  following 
ctiaracter  istics : 

Think  Time  = 300.0  Jiffies 

Processor  Service 

Class  1 = 30.6  jiffies 

Class  2 : 83. b jiffies 

Disk  1/0  Calls 

Class  1 = 10 

Class  2 = 3 


The  results  in  part  1 of  Table  XX  show  that  larger  errors 
are  produced  when  differing  processor  service  time 
r Ovjui rements  are  included.  The  response  time  errors  would 
appear  to  indicate  that  class  1 jobs  are  experiencing  some 
delay  not  anticipated  by  the  model.  This  may  be  correct  or 
some  unexpected  delay  may  be  most  evident  on  the  class  1 
jobs.  Class  1 customers  receive  27.5  jiffies  of  service 
during  the  pass  through  PQ2.  However,  it  was  noted  that  the 
jobs  of  both  classes  were  occassionly  requeued  before 
receiving  a full  quantum  of  service.  The  data  obtained  from 
the  MKrER  program  does  not  provide  the  reason  for 
requeue  ing.  However,  the  most  likely  cause  is  the  expiration 
of  the  job's  in-core  protect  time,  t'lien  this  occurs,  the  job 
is  returned  to  the  end  of  the  queue  even  though  it  may  not 
be  necessary  to  swap  it  out  of  core.  Since  the  jobs  were 
requeued  before  receiving  their  full  quantum,  the  expected 
amount  of  service  required  on  the  last  pass  through  PQ2  was 
caused  to  vary.  Analysis  of  several  runs  of  METHR  data 
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indicated  that  the  jobs  in  this  experiment  received  an  ' 

average  of  2^  jiffies  of  service  on  the  first  pass  through 

PQ2.  This  implies  that  class  1 customers  could  be  expected 

to  make  an  extra  pass  through  the  queue.  The  effect  on  i 

class  2 customers  is  reduced  since  they  require  several  i 

passes  in  any  event.  I 

The  results  of  the  application  of  the  service  time 
adjustment  for  this  experiment  are  shown  in  part  2 of  Table 
XX.  There  is  some  improvement  in  the  errors.  However,  the 
service  time  adjustment  was  based  on  jobs  receiving  an 
allocation  of  30  jiffies  in  PQ2.  If  the  jobs  are  requeued  i 

before  completing  service  because  of  the  expiration  of  in- 
core  protect  time,  then  a new  adjustment  is  needed.  In  order 
to  determine  the  new  adjustment,  it  is  necessary  to  know  the 
amount  of  service  received  in  each  pass  through  the  PQ2  i 

queue.  An  average  value  was  determined  for  this  experiment 
but  the  location  of  the  jobs  within  the  run  queue  when  in- 
core  protect  time  expires  will  vary  with  the  workload. 

The  results  of  experiments  19  and  20  indicate  that 
Closed  Queueing  Network  Models  may  be  applied  with 
reasonable  accuracy  to  the  PECs ystem- 1 0 in  the  single 
processor  configuration.  The  increased  response  time 
highlighted  another  problem  in  applying  this  class  of  model  ^ 

to  complex  computer  systems.  The  expiration  of  the  in-core  ■ 

protect  time  tended  to  have  a randomizing  effect  on  the 
amount  of  service  received  during  the  first  pass  through  the 
Round-Robin  queue  of  PQ2.  Certain  classes  of  customers  wore 


more  sensitive  than  others  to  this  change.  Since  the  amount 
of  time  spent  in  the  run  queue  will  be  partially  dependent 
upon  the  workload,  no  fixed  service  adjustment  is  possible. 
In  order  to  fully  analyze  this  problem  a more  detailed  study 
of  the  scheduling  and  swapping  algorithms  is  needed. 

V.5  Customer  Swapping  Model 

A Closed  Queueing  Network  Model  was  used  to 
represent  the  DECsystem-10  when  customer  swapping  was 
required.  This  model  is  shown  in  Figure  11.  Swapping  is 
required  when  the  total  memory  requirement  of  all  users 
exceeds  the  core  capacity  of  the  system.  Since  the  addition 
of  5 1 2K  words  of  core  memory,  excessive  swapping  has  not 
been  apparent  during  operation  of  the  DFCsystem- 1 0 at  AFAL. 
However,  user  requirements  often  tend  to  expand  to  the 
capacity  of  critical  resources  such  as  memory.  For  this 
reason,  a workload  was  developed  which  induced  swapping  so 
that  the  TRACK  and  METER  programs  could  be  used  to  analyze 
the  effects  on  the  system  performance. 

V . 5 . 1 Ex  pe  t 2 1 

The  workload  for  experiment  21  was  comprised  of  two 

classes  of  10  jobs  each,  with  the  following  characteristics: 

Think  Time  = 300,0  jiffies 

Processor  Service 

Class  1 = 30.6  jiffies 

Class  2 = 83.6  jiffies 

Disk  I/O  Calls 

Class  1 = 10 

Class  2 = 3 
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In  addition  to  the  above  chai-ac  ten  i st  ic  s , each  job  was 
allocated  50K  words  of  core.  The  workload  of  10  jobs  in 


class  1 and  10  jobs  in  class  2 required  a total  of  1000K 
words  of  core.  The  3YSTAT  program  indicated  that  available 
user  memory  could  only  contain  13  of  these  jobs  at  one  time. 
The  remaining  7 would  have  to  be  swapped  onto  a specified 
RPO^  disk  utiit. 

The  additional  load  on  the  Operating  System  imposed 
by  job  swapping  caused  a large  Increase  in  measured 
processor  overhead.  The  overhead  of  CPUI  increased  to  9i, 
while  the  overhead  of  CPUO  was  measured  at  18%.  Ttie  mean 
service  time  requirement  of  each  class  was  calculated  in  the 
same  manner  as  earlier  experiments  without  the  use  of  the 
.service  time  adjustment.  The  average  delay  for  swapping  was 
calculated  from  the  data  of  the  MPTPR  program. 

The  results  of  this  experiment  are  shown  in 
Table  XXI.  Using  the  measured  sw.ipping  delays  proiluced  good 
response  time  errors.  However,  the  availability  of  the 
average  swapping  delay  for  general  workloads  is 
questionable.  Neither  the  TRACK  nor  SY.STAT  programs  provided 
parameters  which  could  be  directly  related  to  the  measured 
values  for  swapping  delay.  The  fact  that  13  customers  could 
occupy  memory  at  one  time  does  not  imply  that  the  Swapper 
Module  does  not  attempt  to  optimise  performance  and  swap  out 
more  jobs  than  than  is  actually  required  for  this 
iJeterm  in  ist  ic  workload.  In  addition,  some  of  the  delay  may 
be  absorbed  as  the  customer  waits  in  the  POI  queue.  This 
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>^ueue  can  tion.ially  be  expected  to  be  small  but  its  length 
will  be  dependent  upon  the  workload. 


Another  factor  to  be  considered  when  swapping  occurs  | 

is  the  interaction  between  user  disk  I/O  operations  and  the 

system  swapping.  A similar  woi'kload  which  required  that  the  i 

synthetic  jobs  write  onto  a RPOU  disk  showed  a drastic  | 

I 

increase  in  measured  I/O  time  when  swapping  was  required.  In 
that  instance,  the  Operating  System  and  the  user  programs 
were  both  using  the  same  disk  control  unit.  The  average 
user  disk  interaction  was  measured  at  1,72  jiffies  when  the 
user  I/O  was  to  a RP03  disk  unit.  When  swapping  was  induced 

and  both  the  Operating  System  and  the  synthetic  jobs  used  , 

the  RPO't  disk  units,  the  average  user  disk  interaction  time 
increased  to  approximately  40  jiffies.  The  size  of  the 

increase  in  I/O  time  warrants  additional  investigation  in  ■* 

later  studies. 

The  representation  of  the  DECsystem- 1 0 shown  in 
Figure  11  has  been  shown  to  produce  relatively  good 
predictions  of  the  system  performance  when  swapping  is 
required.  However,  the  value  of  the  model  is  dependent  upon 
the  availability  of  the  values  for  average  swapping  delay 
for  each  customer. 

I 
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Chapter  VI 

Conclusions  and  Recommendations  for  Future  Study 


VI . 1 Conclusions 


The  Closed  Queueing  Network  Model  provides  a 
flexible  means  of  representing  complex  computer  systems. 
Several  representations  of  the  DECsystem-10  have  been 
presented  in  this  report.  The  analytic  solution  presented  in 
Chapter  III  placed  requirements  on  the  system  being  modeled 
which  are  not  totally  satisfied  by  the  DECsystem-10.  In 
particular,  the  system  does  not  normally  satisfy  the 
limitation  on  service  rates  for  FCFS  service  centers  and  the 
requirement  that  a customer  may  only  occupy  one  resource 
queue  at  a time.  The  validation  of  the  different  models  used 
to  represent  the  system  has  indicated  how  certain  system 
features  and  parameters  of  the  DECsystem-10  affect  the 
accuracy  of  various  model  reijr esentations . 

The  dual  processor  node  model  used  for  experiments  1 
through  6 produced  response  time  errors  ranging  up  to  33%. 
Although  the  accuracy  is  not  as  good  as  would  be  desired  for 
detailed  performance  evaluation  of  the  system,  the 
flexibility  of  Closed  Queueing  Network  Models  permits  many 
configurations  to  be  easily  represented.  The  alternate 
representation  used  for  experiments  7 through  12  indicated 
that  equally  accurate  performance  measures  may  be  obtained 
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with  a simplified  model  using  only  one  Processor  Sharing 
node  to  represent  the  processors.  Finally,  adjustments  were 
made  to  the  service  time  requirements  of  the  customers  in  an 
effort  to  compensate  for  the  limitations  imposed  by  the 
closed  form  solution  of  the  model.  The  results  from  the 
application  of  the  adjusted  service  times,  experiments  13 
through  13,  indicated  that  the  errors  for  worst  case 
workloads  are  improved  while  there  is  a slight  decrease  in 
accuracy  for  other  workloads. 

A Closed  Queueing  Network  Model  was  also  used  to 
represent  the  DECsystem-10  in  a single  processor 
configuration.  The  response  time  errors  were  relatively 
Siiiall  with  the  exception  of  class  2 users  in  experiment  20. 
That  experiment  highliglited  the  influence  of  the  in-core 
protect  time  counter.  Even  though  swapping  was  not  required 
for  the  workload  in  that  experiment,  the  in-core  protect 
time  had  a distinct  effect  on  the  system  performance  for 
specific  customer  classes.  Although  additional  inv est  igat ion 
is  needed  in  this  area,  the  !;iodel  used  for  this 
configuration  provides  a flexible  representation  of  the 
system  and  yields  results  of  reasonable  accuracy. 

Finally,  the  effects  of  swapping  on  system 
performance  were  investigated.  The  results  of  experiment  21 
indicate  that  a delay  node  can  be  used  to  represent  the 
swapping  delay.  However,  the  application  of  the  model 
requires  that  swapping  times  be  taken  from  METER  or 
equivalent  data.  The  collection  of  this  detailevi  information 
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is  generally  not  practical  during  normal  operation.  Kven 
though  the  response  time  errors  were  below  lOt,  the 
requirement  for  the  average  swapping  delay  makes  the 
application  of  this  model  impractical  at  the  present  time. 
Further  study  may  determine  a correlation  between  the 
average  swapping  delay  and  one  or  more  parameters  available 
in  the  TRACK  data.  A relationship  of  this  form  would  permit 
the  application  of  the  model  to  tlie  DKCsystem- 1 0. 

Closed  Queueing  Network  Models  offer  great 
flexibility  for  representing  complex  computer  systems  such 
as  the  DFCs ys tern- 1 0 . The  product  form  solution  to  this 
class  of  model  requires  tliat  assumptions  be  made  concerning 
the  system  to  be  modeled.  This  report  has  shown  that  the 
product  form  solution  for  Closed  Queueing  Network  Models  is 
applicable  to  the  IM-ICsystem- 1 0 if  a high  degi’ee  of  accuracy 
is  not  required. 

VI.  2 Recommend  at- ions  for  Future  .Study 

The  magnitude  of  the  response  time  errors  for  the 
models  presented  in  this  report  arise  from  two  areas; 
limitations  in  the  model  caused  by  assumptions  made  to 
obtain  a closed  forMi  solution  and  limitations  in  the 
construction  of  the  synthetic  workload.  Additional  topics 
for  future  study  include: 

1.  The  limitations  imposed  by  the  assumptions 


required  for  the  model's  solution  need  to  be  further 
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invost  . A rii^orous  niatluMiiatical  invest  igat  ion  should 
be  undertaken  in  order  to  deter:nine  refined  service  time 
adjustment  procedures  and  the  limitations  of  such 
procedures . 

2.  A synthetic  workload  which  more  closely 
approximates  the  assumptions  of  the  model  needs  to  be 
developed  in  order  to  effectively  validate  models  of  this 
nature.  The  SCRIPT  program  permits  the  simultaneous 
execution  of  numerous  programs.  Classes  could  be  modeled  by 
several  synthetic  jobs  loaded  into  the  system  using  several 
"scripts"  for  eaoii  class  of  customer. 

3.  A laore  compi'ohonsi ve  investigation  of  scheduling 
and  swapping  procedures  is  required  to  analyze  the  problem 
of  early  requeueing  of  jobs  in  the  run  queue.  In  addition, 
a sensitivity  study  sliould  be  performed  to  determine  the 
effects  on  performance  resulting  from  changes  in  the 
resource  allocation  parameters. 

Additional  study  of  I/O  features  needs  to  be 
accomplished  to  validate  this  class  of  model  for  the 
complete  system.  This  i n v est i gat i on  should  include  the 
CvTisequences  of  swapping  on  I/O  and  overall  system 
per  formance . 

5.  In  order  for  a model  which  incorporates  customer 
swapping  to  be  of  practical  value,  readily  accessible 
parameters  which  provide  good  correlation  with  observed 
swapping  delays  must  be  identified.  The  identi f icat  ion  of 
these  parameters  will  require  detailed  knowledge  of  the 
Operating  System  .Swapping  routines. 
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Ap  pe  nil  ix  A 

Closed  Queueing  Network  Model  Solution  Technique 


Networks  of  queues  provide  an  flexible  method  of 
modeling  multiprogrammed  and  time-sharing  computer  systems. 
A computer  system  may  be  viewed  as  a network  of  resources 
(terminals,  CPUs,  I/O  devices)  and  a collection  of  customers 
and  their  associated  tasks.  Similar  resources  may  be 
grouped  together  into  a single  class  called  a service 
center.  The  customers  pass  from  one  service  center  to 
another,  receiving  some  amount  of  service  from  each  of  tne 
centers  along  their  route.  If  the  routes  are  such  that  a 
customer  may  neither  enter  not  leave  the  network,  then  the 
network  is  said  to  be  a closed  queueing  network.  In  a time- 
sharing environment,  one  of  the  service  centers  will 
repr '^.ser.t  the  user's  terminals.  The  command  entered  by  the 
user  proceeds  to  the  system  resources  where  it  receives  some 
amount  of  service  and  returns  to  the  terminal  center  where 
it  receives  more  service  in  the  form  of  think  time. 
Following  the  think  time,  another  task  leaves  the  ter::iinal 
service  center  and  proceeds  through  the  network.  The  closed 
(lueueing  network  then  provides  a realistic  model  of  a time- 
sharing computer  system. 

Associated  v;ith  each  service  center  in  the  network 
is  a queue.  Customers  are  assigned  to  a service  center's 
queue  while  they  are  awaiting  or  receiving  service  from  that 


center.  The  state  of  a closeci  queu<‘ing  network  can  then  be 
described  in  terms  of  tlie  number  of  customers  at  each  queue 
in  the  network. 

The  system  chantjes  state  when  a customer  departs  one 
service  center  and  enters  the  queue  of  another.  Jackson 
(Ref  12)  and  Cordon  and  Newell  (REF  13)  have  described 
the  long-run-average  behavior  in  terms  of  an  ecjuilibrium 
state  probability  distribution  function.  At  equilibrium,  the 
rate  of  transition  of  customers  into  a state  from  all  other 
states  must  equal  the  rate  of  transition  of  customers  out  of 
that  state  into  all  other  states.  Tne  equations  describing 
this  feature  are  called  global  balance  equations  by  Chandy 

(Ref  16).  Cordon  and  Newell  used  this  concept  when 

solving  the  balance  equations  for  the  case  of  identical 

customers.  Customers  are  said  to  be  identical  if  their 

routit!g  through  the  network  is  the  same  and  if  their  service 
time  distribution  functions  are  the  same  at  each  service 
center.  The  following  discussion  covers  the  development  of  a 
closed  form  solution  for  tlie  ecjuilibrium  state  probability 
distribution  (Ref  13)* 

While  deriving  the  steady  state  probabilities,  two 
additional  futictiuns  will  be  lielpful.  Let 

0 n . = 0 

e(n.)  ^ ^ (26) 

^ I 1 n ^ > 0 

1 12 


( 


This  function  is  o dimensionless,  multiplicative  factor 
which  will  be  used  to  eliminate  from  the  steady  state 
probabilities  terms  which  cannot  exist.  Thus,  the  final 
result  will  be  simplified  with  the  understanding  that  the 
expression  is  valid  only  for 


n ^ 0 and 


Another  function  wliich  will  simplify  notation  is: 


a . ( n . ) = 
1 1 


n . n ,<  r 

1 i 

r n r 


(27) 


This  function  is  also  a dimensionless,  multiplicative 
factor.  It  will  be  used  to  indicate  the  number  of  Jobs 
receiving  service  at  the  ith  service  center,  with  r 
Identical  servers,  when  there  are  tK  Jobs  in  the  queue. 
Actually  the  factor  may  be  any  positive  function,  but  the 
use  of  this  particular  will  clarify  its  use  in  the 

development  of  the  steady  state  probi'aility  distribution. 

The  network  will  be  characterized  by  the  following 
parameters : 

N = the  number  of  customers  (tasks  or  Jobs) 
receiving  service  in  the  network. 

p,  j = the  probability  tliat  a customer  leavirq;  the 
ith  service  center  will  proceed  to  the  Jth  service  center. 


113 


T 


M = the  number  of  service  centers  in  the  network, 
r^  = the  number  of  identical  servers  at  the  ith 
service  center. 

= the  departure  rate  from  a server  in  the  ith 

service  center.  Expressed  in  jobs/sec. 

A.  = the  arrival  rate  at  the  ith  service  center.  As 
1 — 

with  the  departure  rate,  expressed  in  jobs/sec. 

P(  n ^ ,n  2 I . . . , = the  joint  steady  state  probability 


that  there  are  n^  customers  at  the  ith  service  center.  These 

probabilities  should  only  be  defined  for 

Consider  the  rate  of  transition  out  of  a state 

( n , n - : 

I d M 


E 


k=1 


6 (n^)a|^(n,^)  Mj^PCn^  ’*^2’  ‘ ' ' 


(28) 


This  formula  indicates  that  customers  exit  from  state 
( n pn^,  • . -nji)  throush  the  kth  center  as  long  as  there  is  at 
least  one  customer  present  at  that  center.  If  there  are  no 
customers  at  the  kth  center,  no  departures  are  possible  and 
the  function  will  set  the  corresponding  term  equal 

to  zero  for  the  summation. 

The  rate  of  transition  into  a state  ( n , , n _ , . . . n,.) 

I d M 


will  be  expressed  as: 


P(n 


,n  . + 1 , 


i = 1 


k=1 

k^i 


E 

i = 1 


E ( n . ) a . ( n . ) 
111 


M.  p . . 

Ill 


P(n 


1’“^' 


(29) 


This  expression  indicates  that  a transition  into  a state 
( n ^ ,n  2 I • • occurs  from  state  ( n ^ , . . . . n 1 , . , n 1 , . . n 1^^) 

whenever  a customer  completes  service  at  the  ith  center  and 
then  proceeds  to  the  kth  center.  then  n = 0,  the  state 

K 

( n j . n 1 , . . . n 1 , . . . n cannot  exist.  However,  for 

mathematical  convenience,  P ( n ^ . n 1 , . . n 1 , . . n will  be 
allowed  to  exist  and  E(n,)  will  eliminate  that  term  from 
the  summation  (Ref  22; 2^0) 

The  (lueue  dependent  function  a (n  ) must  be  handled 

K K 

differently  when  i = k.  In  this  case,  a customer  departs  the 

ith  service  center  and  immediately  returns  to  that  center. 

Therefore,  the  queue  dependent  departure  rate  must  be  based 

on  n.  customers  rather  tlian  n.  + l customers. 

1 1 

At  equilibrium,  the  rate  of  transition  out  of  any 
state  is  equal  to  the  rate  of  transition  into  that  state. 
Therefore,  the  steady  state  probabilities  must  satisfy: 


Z •'"> 


M,^P(n  j .n^ 


M M 


^ S 6 (n^)a.(n  .+  1 ) M.p 


i'^ik  ' • • '"k“^  ’ 


,tli+1  I • • 


ri 

^ 6(n.)a,(n.)  M.p.  . i>(  n ^ , ri^,, . . . 


Now  consider  the  equilibrium  condition  for  a particular 


service  center: 


Number  of  Jobs  Departing,  = Number  of  Jobs  Arriving 


Note  that  jobs  may  only  (iepart  from  a service  center  if  the 


center  is  being  utilized  and  that  it  may  be  utilized  only 


when  at  least  one  job  is  present  in  its  queue.  The 


expression  for  the  loft  side  of  the  equation  is 


* (Active  Time),  while  the  right  hand  side  is 

* (Total  Time).  Reordering  tei'ins  yields  the  definition 


for  resource  utilization: 


Percent  Utilization  = (Active  Time) /(Total  Time)  = A / n, 

1 i 


M 


1 


r 


( 


L 


This  ratio  is  also  called  the  traffic  intensity  and  is 
sometimes  given  units  of  erlangs  in  honor  of  A.  K.  Erlang,  a 
pioneer  in  congestion  theory  (Ref  23;'40).  The  are 

parameters  of  the  queueing  network  and  will  remain  constant. 
The  will  vary  with  the  workload  N.  It  is  assumed  that 

changes  in  the  number  of  jobs  in  the  network  will  not  affect 
the  relative  utilizations  of  the  service  centers.  This  is 
intuitively  acceptable  when  all  the  customers  are  identical. 
However  when  customers  have  differing  service  requirements 
at  a First  Come  First  Serve  queue,  then  the  assumption  is  no 
longer  valid;  as  is  discussed  in  Chapter  III. 

For  centers  with  identical  servers,  the  utilization 

is : 


X .. 


K. /u  . 


a.(n.)  a.(n.) 


(31) 


Where  is  the  utilization  of  an  equivalent  single  server 
center.  Doth  X and  are  functions  of  the  workload,  N. 
Wlien  the  service  times  are  state  dependent  and  the  function 
a^(n^)  is  an  arbitrary  positive  function,  the  utilization  of 
tlie  servers  at  the  ith  center  is  no  longer  expressed  by 
Eq  (31).  In  this  instance,  the  concept  of  an  individual 
server  within  the  service  center  is  no  longer  clearly 
defined . 

In  order  to  solve  the  balance  equations,  the  method 
of  seperation  of  variables  will  be  used.  As  stated  earlier, 
it  is  assumed  that  the  utilizations  of  the  service  centers 


1 I? 


are  functions  of  the  network  parameters  p . and  m , and  of 

ij  i 

the  number  of  customers  in  the  network.  Once  an  equilibrium 
condition  is  reached,  the  remain  constant.  These 

utilizations  are  then  used  to  define  a new  function  as 
follows : 

. M n . 

Q(n  ,n^,  . . . ,n  ) = J_  p-I  (X  .)  ^ (32) 

' " C(N)  1 1 ^ 

i=  1 


khere  C is  a normalization  constant  which  will  later  be  used 
to  adjust  for  the  dependence  of  the  X^  upon  the  number  of 
customers  in  the  network. 

For  notational  convenience  also  define: 


11,(0) 


I3,(n) 


= 1 


ak(n)  Bj^(n-I) 


k= 1 ,2. . .M 
n=1 ,2.  . .N 


(33) 


These  functions  permit  the  following,  chanije  of  variables  to 
be  made  in  the  balance  equations: 


Q ( n ^ , n n ^^) 

[i 

k = I 


(3n) 
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Note  that 


P(n. 


.n^  + 1 


•n.,) 


Q(n^  , . .n|^-1  . .n^  + 1 . .n|;^) 

( n^  + 1 ) 


M 

II 


Qj(nj) 


(35) 


Substituting  into  the  balance  equations  and  canceling  the 
terms  yields  the  following  results: 


Z «("k 


k=  1 


*'^2’  • • 


M M 


= 22  S(nj^)a^(n^)  M.  p^.  Q(  n ^ . n^- I , . . n 1 , . . n,,, ) 


(36) 


i=  1 k=  1 


substituting  for  Q(-)  yields: 


M ^ M 

2 S("k)^k("k)"k 

k=  1 


n 


i=  1 


M M 


^ M 


= Z Z E(.v)n^(npM;P(,  (Xj/X,)  rixjj 
i=lk=1  Lj=) 


(37) 


which  reduces  to: 


k:  1 


M M 

1= 1k=  1 


€(n^)a^(n^)Pk  ' 


0 


(38) 


TOU'  liUl".  J 


-J  MEL"*' 


\ 
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Since  all  the  jobs  could  bo  at  any  of  the  service  centers  at 
one  time,  the  following  must  hold: 

M 

" ‘39) 

1=  1 

From  the  earlier  definition  of  X.  , this  can  be  written  as: 

1 

M 

X A p = A k=  1 ,2 M (J)0) 

i = 1 i ik  k 

Tliis  is  not  an  independent  set  of  equations.  The 
can  only  be  determined  within  a multiplicative  constant.  The 
interpretation  of  the  set  of  equations  is  that  the  arrival 
rate  to  the  kth  center  is  equal  to  the  weighted  sum  of  the 
arrival  rates  of  the  service  centers  which  could  possibly 
send  it  a customer. 

If  the  network  is  modeling  a time-sharing  computer 
system,  then  the  relative  arrival  rates  will  be  of  intei'est. 
Let  the  user  terminals  be  represented  as  workcenter  number 
one.  Then  a^  is  the  relative  arrival  rate  at  service 

center  i with  respect  to  the  user  terminals.  This  indicates 
that  a customer  will  visit  the  ith  service  center  A 

i 1 

times  between  the  time  it  loaves  the  terminal  until  it 
returns. 

Define  o ^ as  the  expected  number  of  visits  to  the 
ith  center  per  interaction  with  the  user  terminal  service 
center.  This  permits  the  last  evjuation  to  be  written  as: 
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p 


i=1 


e.  p. , 
Ilk 


"k 


k=1 ,2, . . .M 


(41) 


Now  the  relative  number  of  visits  to  each  center  can  be 
computed  as  well  as  the  relative  utilizations  for  identical 
servers,  . Rewriting  the  equilibrium  state  probability 
distribution  function  yields: 


1 

P ( n , n , . . . n ) = 

12  M c(N) 


n 


i=  1 


1 1 


(42) 


i ) 


Since  the  X^  terms  are  utilizations,  and  the  ( n^  ) terms 
are  products  of  multiplicative  factors,  there  are  no  units; 
as  expected  for  a probability  distribution.  Note  that: 

n ^ 

n i(n  i)  = n .(k)  ('13) 

k=  1 ^ 

and  let 

M.(k)  = M.a  .(k)  (44) 

J J J 

khere  M.(k)  is  the  departure  rate  from  the  j^h  center  when 
d 

there  are  k customers  at  the  center.  Tlien: 


P ( n J , n ^,  . . . n i^j) 


'UN) 


F,(n,)P'  (n.,)  . . .F,.(n,.) 


r r‘  2'  2 


M M 


(4‘3) 
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where 


F. (n. ) = (46) 

1 1 

Mi(  j) 

The  F tj  ) terms  are  unnormalized  probabilities  that  service 
center  i has  n^  customers  in  its  queue. 

The  normalization  constant  is  now  written  as  C(N)  to 
indicate  its  dependence  upon  the  workload,  N.  C(N)  must 
normalize  the  product  on  the  right  hand  side  of  the  equation 
to  insure  that  the  probabilities  sum  to  one  over  all 
possible  states.  Therefore: 


-1 


C(N) 


Z 


F ^(n^)F2(n2)  • • ‘F.iCnj^^) 


(47) 


all  states 
n +n_-t- . . . +n,,=  N 

1 d rl 
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